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^described.  This  report  is  the  fourth  in  series  Technical  Repcrt  H-T5-** 
published  under  the  general  title  "Los  Angeles  and  Long  Beach  Harbors  Model 
Study.’  A model  distortion  ratio  of  1:1<  and  a vertical  scale  ratio  of 
1:100  (l:L00  horizontal  scale  ratio)  vere  selected  to  minimize  the  model 
area  and  to  provide  a vertical  scale  ratio  where  accurate  model  wave  mea- 
surements could  be  assured.  The  model  reproduces  approximately  253  square 
miles  of  prototype  area.  Model  limits  were  selected  to  enclose  a strong 
convergence  zone  seaward  of  the  harbor  breakwater.  Adjustment  of  the 
initial  generated  wave-front  position  in  the  model  was  necessary  to  com- 
pensate for  model  distortion  and  the  maximinn  model  depth  of  -300  ft  mllw. 

rhe  electrohydraulic  wave  generator  used  in  the  model  study  is  com- 
posed of  iL  separate  15-ft  units  for  a total  wave  generator  length  of 
210  ft  and  is  capable  of  (a)  generating  waves  with  a prototype  period 
ranging  from  15  to  600  sec,  (b)  generating  a wave  with  small  variation  in 
period  and  height  between  units,  (c)  defining  resonant  response  over  a 
narrow  period  band,  (d)  generating  a variable  wave  height  along  the  wave 
front,  and  (e)  generating  a curved  wave  front Performance  requirements 
for  the  wave  generator  include:  (a)  the  ampliwde  error  of  the  wave-paddle 

stroke  for  each  unit  must  be  less  than  1 percent  and  (b)  the  phase  lag  of 
the  wave-paddle  response  from  the  electronic  command  signal  between  each  of 
the  iL  units  must  be  less  than  1*  degrees.  \ 

The  four  subsystems  of  the  automated  model  qpta  acquisition  and  con- 
trol system  include:  (a)  digital  lata  recording  and  controls,  (b)  analog 

recorders  and  channel  selection  circuits,  (c)  wave  sensors  and  interfacing 
equipment,  and  (d)  wave  generator  units  and  control  equipment. 

Model  wave  record  analyses  include  the  capacity  for  determination  of: 
(a)  average,  root  mean  square,  and  significant  wave  height  and  (b)  average 
period.  Procedures  for  Fourier  analysis  and  least-squares  harmonic  analysis 
of  the  wave  record  are  also  available. 

The  hydraulic  model  can  be  used  in  a predictive  mode  to  evaluate 
various  proposed  modifications  to  optimize  the  modifications  from  both  a 
cost  and  functional  standpoint.  Potential  future  uses  of  the  model  in- 
clude studies  of  (a)  dispersion  of  effluent  from  specific  outfalls  in  the 
harbor  complex,  (b)  tsunami  waves,  (c)  effects  of  various  coastal  struc- 
tures from  Point  Fermin  to  Huntington  Beach,  (d)  effects  of  artificial  off- 
shore islands  or  landfill  seaward  of  the  existing  harbor  breakwater  on 
tidal  circulation  and  wave  conditions,  (e)  tidal  circulation  in  areas  such 
as  Huntington  Beach,  Bolsa  Chica  Wildlife  Reserve,  Alamitos  Bay,  etc., 

(f)  effect  of  various  future  harbor  modifications  on  tidal  circulation  and 
wave  conditions,  (g)  studies  of  the  physical  environment  in  the  harbor  and 
surrounding  areas,  and  (h)  wave  and  tidal  circulation  studies  of  the  pro- 
posed Downtown  and  Queensway  Marinas. 
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PREFACE 


The  Los  Angeles  and  Long  Beach  Port  Authorities  plan  to  construct 
additional  harbor  basins  and  dredge  deeper  channels  and  harbor  ^reas  to 
meet  growing  future  demands  for  expansion  of  ship  mooring  facilities. 

The  Los  Angeles  and  Long  Beach  Harbors  model  study  was  undertaken  to 
investigate  (a)  incidence  of  harbor  oscillations  adverse  to  ship  moor- 
ing, (b)  tidal  circulation  characteristics  of  the  existing  harbor  and 
proposed  expansions,  (c)  optimum  plan  for  future  expansions,  and  (d)  ef- 
fects of  future  expansions  on  existing  harbor  conditions.  This  report, 
the  fourth  in  a series  to  be  published  under  the  general  title  "Los 
Angeles  and  Long  Beach  Harbors  Model  Study,"  considers  design  criteria, 
wave  generator  design  and  control,  and  data  acquisition  for  the  hydrau- 
lic model. 

Project  administration  and  funding  were  provided  by  the  U.  S. 

Army  Engineer  District,  Los  Angeles  (SPL),  under  project  management  of 
Messrs.  J.  Chapman,  H.  Converse,  and  T.  Nizinski  and  under  the  general 
direction  of  Messrs.  G.  Fuquay,  former  Chief  of  the  Engineering  Divi- 
sion, T.  Nishihara,  Chief  of  the  Engineering  Division,  and  C.  H.  Fisher, 
Chief  of  the  Coastal  Resources  Branch.  COL  R.  J.  Malley,  CE,  COL  J.  V. 
Foley,  CE,  and  COL  H.  G.  Robinson,  CE,  were  District  Engineers  of  SPL 
during  the  course  of  this  study.  General  project  administration  for 
the  U.  S.  Army  Engineer  Division,  South  Pacific,  was  provided  by 
Messrs.  0.  F.  Weymouth,  0.  T.  Magoon,  J.  W.  Gerhart,  and  A.  E.  Wanket. 

The  data  analyses  were  conducted  in  the  Hydraulics  Laboratory  of 
the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  under  the 
general  supervision  of  Messrs.  H.  B.  Simmons  and  F.  A.  Herrmann,  Chief 
and  Assistant  Chief,  respectively,  of  the  Hydraulics  Laboratory.  This 
report  was  prepared  by  Mr.  D.  G.  Outlaw,  Harbor  Wave  Action  Branch 
(HWAB),  Dr.  D.  L.  Durham,  HWAB,  Mr.  C.  E.  Chatham,  Chief  of  HWAB,  and 
Dr.  R.  W.  Whalin,  Chief  of  the  Wave  Dynamics  Division. 

Model  construction  was  supervised  by  Mr.  C.  W.  Brasfeild  with  the 
assistance  of  Mr.  R.  R.  Bottin.  Two-dimensional  breakwater  and  harbor 
side-slope  wave  tests  were  conducted  by  Messrs.  Chatham  and  L.  A.  Barnes, 
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Civil  Engineering  Technician.  The  prototype  wave  refraction  analysis 
was  conducted  by  Mr.  Chatham  and  the  model  refraction  analysis  by 
Mr.  Outlaw. 

Directors  of  WES  during  the  model  design  and  the  preparation  and 
publication  of  this  report  were  BG  E.  D.  Peixotto,  CE,  COL  G.  H.  Hilt, 
CE,  and  COL  John  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown 
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CONVERSION  FACTORS,  U,  S.  CUSTOMARY  TO  METRIC  (Sl) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Sl)  units  as  follows: 


Multiply 


To  Obtain 


inches 

feet 

miles  (U.  S.  statute) 

miles  (U.  S.  nautical) 

fathoms 

square  feet 

square  miles  (U.  S. 
statute) 

acres 

cubic  yards 

gallons  per  day 

pounds  (force) 

feet  per  second 

feet  per  second  per  second 

degrees  (angle) 

Fahrenheit  degrees 


25. h 
0.301*8 

1.60931+1* 

1.852 

1.8288 

O.O929030I* 

2.589988 


millimetres 

metres 

kilometres 

kilometres 

metres 

square  metres 
square  kilometres 


1*01*6.856 

0.761*551+9 

3.T85I+I2 
1*.  1*1*8222 
0.301*8 
O.30I+8 
0.0171*5329 
5/9 


square  metres 
cubic  metres 

cubic  decimetres  per  day 
newtons 

metres  per  second 

metres  per  second  per  second 

radians 

Celsius  degrees  or  Kelvins* 


• To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (f)  read- 
ings, use  the  following  formula:  C = (5/9)(F  - 32).  To  obtain 

Kelvin  (k)  readings,  use:  K = (5/9)(F  - 32)  + 273.15. 
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LOS  ANGELES  AND  LONG  BEACH  HARBORS  MODEL  STUDY 


MODEL  DESIGN 
PART  I : INTRODUCTION 

Background  and  Model  Study  Ob.lectlves 

1.  The  ports  of  Los  Angeles  and  Long  Beach  are  located  on  San 
Pedro  Bay  along  the  southern  coast  of  California.  A location  map  and 
the  harbor  configuration  are  shown  in  Figure  1.  Development  of  the 
harbors  and  historical  surge  activity  are  reviewed  in  detail  as  a 
portion  of  a study  completed  by  Science  Engineering  Associates  (SEA)^ 
for  the  U.  S.  Army  Engineer  District,  Los  Angeles  (SPL),  in  July  I968. 
The  four  major  objectives  of  the  Los  Angeles  and  Long  Beach  Harbors 
study  conducted  by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  are  to: 

Determine  incidence  and  severity  of  troublesome  oscilla- 
tions in  the  present  harbor  complex. 

b.  Investigate  tidal  circulation  characteristics  of  the 
present  and  proposed  harbors. 

c_.  Determine  the  optimiom  plan  for  future  expansions  to  pro- 
vide safe  and  economical  berthing  areas. 

d.  Analyze  the  effect  proposed  expansions  will  have  on 
existing  harbors. 

The  study  results  are  being  published  in  a series  of  reports  under  the 
general  title  "Los  Angeles  and  Long  Beach  Harbors  Model  Study."  This 
report  is  the  fourth  in  the  series  and  covers  model  scale  and  size, 
wave  generator  design,  and  model  data  acquisition  and  analysis  proce- 
dures, primarily  for  wave-induced  harbor  oscillations.  Model  operation 

and  appurtenant  equipment  for  tidal  circulation  studies  are  discussed 
2 

in  Report  5* 

2.  To  achieve  the  first  objective  and  acquire  prototype  tidal 
data  for  the  second  objective,  the  WES  approach  was  to: 

a.  Acquire  prototype  wave  data  over  a 1-yr  period. 

b.  Record  observations  of  ship  motion. 

£.  Catalog  ships  using  Los  Angeles  and  Long  Beach  Harbors. 
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d.  Conduct  analytical  investigations  of  moored  ship  response. 
Perform  extensive  analyses  of  prototype  wave  data. 

_f.  Attempt  to  correlate  ship  motion  with  wave  height  and 
frequency. 

£.  Collect  prototype  tidal  height  and  current  data  to  verify 
model  tests  of  tidal  circulation  within  the  harbor  complex. 

The  prototype  data  collection  program  for  waves,  ship  motion,  and  tides 

3 

IS  described  in  Report  1 and  results  of  the  wave  data  analyses  are 

described  in  Report  3 of  this  series.  Ship  motion  data  were  published 

5 2 

in  Report  2 and  tidal  data  in  Report  5- 

3.  The  purposes  of  the  physical  model  study  were  to  investigate 
harbor  oscillation  characteristics  and  tidal  circulation  for  the  exist- 
ing harbor  and  to  optimize  design  of  the  proposed  harbor  modifications 
to  provide  safe  and  economical  berthing  conditions  while  maintaining  or 
improving  circulation  in  the  harbors. 

Several  previous  model  studies  of  San  Pedro  Bay  have  been 
conducted.  The  first  model, ^ built  in  1927,  was  of  a relatively  small 
section  of  Los  Angeles  Harbor  and  was  used  primarily  to  study  the  ef- 
fect of  various  proposed  modifications  on  surge  conditions  in  East 

1 6 T 0 9 

Channel.  Other  model  studies  ’ ’ ’ ’ at  VTES  and  at  the  California 
Institute  of  Technology  were  conducted  to  study  proposed  plans  for  con- 
struction of  a model  for  the  Naval  Base  on  Terminal  Island,  a proposed 
naval  supply  depot  to  be  located  near  Point  Fermin,  and  the  Southeast 
Basin  in  Long  Beach  Harbor.  In  recent  years,  plans  have  been  developed 
for  additional  expansion  of  the  harbors  (described  in  paragraphs  13-18) 
that  would  significantly  modify  the  shape  of  the  cuter  harbors  and 
increase  the  depth  in  some  areas  to  accommodate  deep-draft  vessels.  As 
indicated  in  Reference  2,  proposed  plans  are  so  extensive  that  a careful 
examination  of  their  effects  on  existing  harbor  and  expanded  harbor 
facilities  is  required  to  ensure  tViat  the  optimal  cost  effective  ex- 
pansion plan  is  developed  in  order  to  minimize  the  potential  for  un- 
desirable effects  which  could  prove  either  irreversible  or  extremely 
costly  to  correct.  Wave-induced  harbor  oscillations  can  be  of  suffi- 
cient magnitude  to  create  navigation,  berthing,  and  mooring  problems 


T 


due  to  their  effect  on  ship  motion.  Should  existing  tidal  circulation 
be  impaired,  the  potential  exists  for  insufficient  pollutant  dissemina- 
tion and  undesirable  water  quality.  A physical  hydraulic  model  is  an 
excellent  method  for  quantitatively  investigating  harbor  oscillations 
and  tidal  circulation.  Advantages  of  the  physical  model  are  that  hy- 
drodynamic phenomena  can  be  scaled  directly  from  model  to  prototype^ 
(provided  it  is  designed  properly)  and  modifications  to  proposed  harbor 
plans  can  readily  be  evaluated.  Improper  model  design  can  result  in 
the  introduction^  of  spurious  boundary  effects  and  unnatural  wave 
heights  and  mode  shapes  in  the  model. 

5.  In  the  SEA  study,  a physical  model  covering  an  area  of  ap- 
proximately 9-  by  l6-nautical  miles*  with  a l:U-scale  distortion  was 
recommended  for  wave  periods  greater  than  2 min  and  less  than  1 hr.  A 
second  model  of  approximately  6-  by  11-nautical  miles  with  a l:2-scale 
distortion  at  a larger  scale  was  suggested  for  wave  periods  between 

15  sec  and  2 min.  In  the  second  model,  the  effect  of  model  distortion 
on  the  time  scale  ratio  at  lower  wave  periods  would  be  reduced.  General 
model  laws  and  scales  for  modal  shapes  and  frequencies  of  oscillations, 
refraction,  shoaling,  and  diffraction  also  were  discussed  in  the  SEA 
study  for  distorted  physical  models. 

6.  After  review  of  the  preceding  physical  model  recommendations, 
the  feasibility  of  using  one  model  with  a larger  horizontal  scale  was 
considered  at  WES  and  a single  model  (Figure  2)  was  found  to  be  feasi- 
ble for  investigation  of  both  harbor  oscillation  phenomena  and  tidal 
circulation.  In  Part  of  this  report,  the  effects  o.  wave  refraction, 
diffraction,  viscous  ix-iction,  wave  reflection,  wave  transmission 
through  the  harbor  breakwaters,  wave  filters,  and  wave  absorbers  are 
considered  in  order  to  determine  model  scales  and  limits.  In  the  re- 
maining sections  of  the  report,  wave  generator  design  and  model  data 
acquisition  and  analyses  are  described. 


• A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  is  presented  on  page  U. 
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Figure  2.  Los  Angeles-Long  Beach  Harbors  model 


Description  of  Los  Angeles  and  Long  Beach  Harbors 

7.  Los  Angeles  and  Long  Beach  Harbors  occupy  a major  part  of 
San  Pedro  Bay  and  are  separated  by  the  city  boundaries  as  shown  in 
Figure  3.  The  harbor  is  370  nautical  miles  southeast  of  San  Francisco 
Bay  and  is  95  nautical  miles  northwest  of  San  Diego  Bay.  San  Pedro  Bay 
is  protected  on  the  west  by  the  Palos  Verdes  San  Pedro  Hills  and  Point 
Fermin  but  is  exposed  on  the  south  and  southeast.  Santa  Catalina  Is- 
land and  San  Clemente  Island  partially  shelter  the  harbor  from  the 
southwest.  The  offshore  topography  of  the  ocean  bottom  near  San  Pedro 
Bay  is  shown  in  Figure  As  indicated  on  the  figure,  the  continental 
shelf  extends  outward  in  a convex  shape.  Near  the  100-fathom  contour 
line,  the  shelf  drops  off  sharply  toward  the  San  Pedro  Channel,  reaching 
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Figure  3.  Location  of  city  boundary  and  various  channels  and  basins 
in  the  Los  Angeles-Long  Beach  Harbors 


Figure  1* . Underwater  contours 
seaward  of  Los  Angeles-Long 
Beach  Harbors  (contours  in 
fathoms ) 


depths  in  excess  of  500  fathoms.  A submarine  canyon  also  intersects  the 
continental  shelf  on  each  side  of  the  bay,  and  the  canyons  extend 
within  the  model  limits  as  shown  in  Figure  U. 

8.  The  harbor  is  protected  on  the  seaward  side  by  an  8-mile-long 

stone  breakwater  extending  from  Point  Fermin  eastward  past  the  Long 
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Beach  Harbor  (Figure  3).  The  breakwater  consists  of  three  sections 
originally  built  approximately  along  the  50-ft  contour  line.  The  first 
to  be  constructed  was  the  San  Pedro  breakwater  from  Point  Fermin  to 
Angels  Gate,  the  Los  Angeles  Harbor  entrance  channel;  the  second  was 
the  Middle  breakwater  extending  from  Angels  Gate  to  Queens  Gate,  en- 
trance channel  to  Long  Beach  Harbor;  and  the  third  was  the  Long  Beach 
breakwater. 

9.  The  11 , 15?-ft-long  San  Pedro  breakwater  is  of  stone-mound 
•nstructi  n with  a cap  course  of  granite  blocks  to  an  elevation  of 
14  f*  ab<)ve  m<=*an  lower  low  water  (mllw).  The  Middle  and  Long  Beach 
br<'<ii'wit  (.rs  -tlst  are  of  stone-mound  construction  but  have  a sand  core 
''a{p-:  wi'h  a relatively  impervious  clay;  the  crown  elevation  of  these 
'w  treaKwat ers  also  is  +l4  ft  mllw.  The  impervious  core  in  the  Middle 
treaKWH’er  expends  approximately  24  ft  above  the  harbor  bottom  (-26  ft 
m..w  ar.  1 the  core  of  the  Long  Beach  breakwater  is  26  ft  above  the 

bi  ’’  m ,-.'4  ft  mllw).  Typical  cross  sections  through  the  three  break- 
wa’ers  are  showi.  in  Figure  5. 

. jepths'"^  in  Los  Angeles  Harbor  vary  from  l6  to  50  ft.  Con- 
’ ..iiig  iepths  for  various  channels  and  basins  are: 

Channel/Basin Controlling  Depth,  ft  mllw 


Fr. trance  channel 

4o 

Milsectiivn  of  entrance  channel 

51* 

I'urning  basin  (outer  harbor) 

51* 

Los  Angeles  channel 

35 

Turning  basin  (inner  harbor) 

35 

Fast  Basin  channel  and  Cerritos  Channel 

35 

West  Basin 

35 

East  Basin 

35 

• Federal  project  depth  of  4o  ft;  dredged  to  51  ft  by  Port  of  Los 
Angeles. 
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MIDDLE  BREAKWATER 


LONG  BEACH  BREAKWATER 
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io 
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Figure  5.  TVplcal  cross  sections  of  the  Los  Angeles-Long 
Beach  Harbors  breakwater 


In  Long  Beach  Harbor,  depths  vary  from  l8  to  TO  ft.  Controlling  depths 
for  various  channels  and  basins  are: 


Channel /Basin 

Controlling  Depth,  ft  mllw 

Entrance  channel 

62 

Cerritos  Channel 

50 

Turning  basin 

55 

The  Federal  project  depth  for  the  Long  Beach  turning  basin  and  channels 
is  35  ft  mllv,  but  the  depths  have  been  increased  by  the  Port  of  Long 
Beach  to  accommodate  deep-draft  vessels.  Considerable  subsidence  also 
has  occurred  in  years  past  due  to  extraction  of  oil  from  the  area. 

11.  Tidal  circulation,  freshwater  discharge,  salinity,  and 

thermal  stratification  conditions  in  the  existing  harbors  are  discussed 

in  Reference  2.  A brief  suitmary  of  the  data  from  Reference  2 pertinent 

to  model  design  is  as  follows: 

Tides  in  San  Pedro  Bay  are  of  a mixed  type  (two  unequal 
tides  per  day)  with  a mean  tidal  range  of  3.8  ft  and  a 
diurnal  (mean  higher  high  to  mean  lower  low)  range  of 
5.^  ft.  Tidal  datum  is  mllw  which  is  2.8  ft  below  mean 
sea  level. 

b.  Tidal  currents  in  the  bay  are  rather  weak  with  a normal 
maximum  current  velocity  of  approximately  1 fps.  Wind- 
generated currents  can  be  of  the  same  order  of  magnitude. 

c_.  Freshwater  discharges  are  limited  to  intermittent  storm 
runoff  primarily  through  the  Dominguez  Channel  and  the 
Los  Angeles  River  and  a few  other  effluent  discharges. 

d.  Salinity  of  the  harbor  is  essentially  uniform  and  very 
close  to  that  of  the  surrounding  coastal  waters,  approxi- 
mately 33  to  3^  ppt  total  salts.  A surface  layer  of  low 
salinity  may  form  in  individual  basins  due  to  storm  run- 
off from  rains,  but  these  conditions  are  relatively  rare 
and  do  not  persist. 

£.  Low-salinity  effluents  include  cannery  wastes  in  the 

vicinity  of  Fish  Harbor  and  discharge  from  the  Terminal 
Island  Sewage  Treatment  Plant.  The  effect  of  numerous 
other  small  effluent  discharges  in  the  harbor  is  not 
significant. 

£.  Cooling  water  discharges  in  Los  Angeles  Harbor  are  prin- 
cipally due  to  the  Union  Oil  refinery  which  averages 
about  26  mgd  with  a maximum  flow  of  about  30  mgd,  and 
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the  City  of  Los  Angeles'  harbor  steam  plant  located 
adjacent  to  West  Basin.  The  plant  withdraws  water  from 
Los  Angeles  Harbor  slip  5 and  discharges  into  the  north- 
east corner  of  West  Basin  with  an  average  flow  rate  of 
78  mgd.  The  plant  operates  at  intervals  dictated  by 
peak  electrical  demand  and  has  a capacity  flow  rate  of 
397  mgd  with  a 12-15°F  temperature  rise. 

£^.  Thermal  stratification  in  San  Pedro  Bay  ranges  from  mild 
seasonal  temperature  gradients  to  strong  local  gradients 
near  cooling  water  discharges.  Ambient  surface  tempera- 
tures average  about  50°F  during  the  winter  and  68°F 
during  the  summer. 

An  additional  cooling  water  discharge  is  anticipated  from  a return  to 
partial  service  of  the  Southern  California  Edison  Company  Generating 
Station  on  Terminal  Island  with  a maximum  flow  rate  of  about  80  mgd. 

A second  standby  plant  may  increase  the  maximum  flow  rate  by  about 
123  mgd  at  times.  The  cooling  water  intake  and  discharge  are  in  the 
Long  Beach  back  channel. 

2 

12.  The  average  depth  of  the  entire  harbor  is  approximately 
36  ft  below  mllw;  neglecting  the  inner  harbor  area,  the  average  depth 
is  approximately  38  to  39  ft.  Navigation  channels  and  ship  berthing 
areas  usually  exceed  the  average  depth,  as  discussed  in  paragraph  10. 
Average  depths  for  various  sections  of  the  harbor  are  given  in 
Reference  2. 


Los  Angeles  and  Long  Beach  Harbors 
Proposed  Development  Plans 

13.  The  proposed  Federal  project^^  includes  dredging  ship  chan- 
nels and  turning  basins  in  Los  Angeles  Harbor  and  deposition  of  the 
dredged  material  in  parts  of  the  harbor  where  new  land  areas  are  planned. 
The  proposed  Federal  project  will  require  dredging  of  approximately 
10,000,000  cu  yd  and  will  provide  about  I87  acres  of  landfill.  The 
dredged  channel,  shown  in  Figure  6,  is  needed  in  the  Los  Angeles  Harbor 
to  permit  use  of  the  port  by  larger  general  cargo  bulk  and  container 
vessels.  Diking  for  the  landfill  would  be  provided  by  local  interests. 


ll* 
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Figure  6,  Dredging  areas  in  the  Los  Angeles-Long  Beach  Harbors  for  the  proposed  Federal  project 


lU.  Recommendecl  dimensions  of  waterways  in  the  Los  Angeles  Harbor 
Federal  project  are: 

Depth , 


Width,  ft 

Length,  ft 

ft  mllw 

Entrance  channel 

1000 

5,500 

Los  Angeles  channel 

750 

12,500 

I45 

Turning  basin  (inner  harbor) 

1350 

1,650 

East  Basin  channel 

350 

6,000 

»45 

West  Basin 

350-1350 

’,800 

East  Basin 

350-1350 

vara  able 

i 

15*  As  indicated  in  Reference  11,  the  propoL  Federal  project 
is  concerned  only  with  meeting  the  immediate  needs  for  improvement  of 
Los  Angeles  Harbor.  Further  investigation  of  the  Los  Angeles-Long  Beach 
Harbors  area  is  authorized  to  determine  if  there  is  a need  for,  and  a 
Federal  interest  in,  future  navigation  improvements. 

16.  More  extensive  proposed  harbor  development  plans  are  being 
considered  by  the  Port  of  Los  Angeles  and  the  Port  of  Long  Beach  to 
provide  for  supertanker  facilities  and,  for  the  Port  of  Los  Angeles, 
liquified  natural  gas  (LNG)  facilities.  Exact  configurations  have  not 
been  determined  by  the  ports,  and  a number  of  different  configurations 
and  construction  phases  are  to  be  tested  in  the  hydraulic  model. 

17.  The  planned  LNG  facility  in  Los  Angeles  Harbor  will  be  con- 
structed on  a landfill  area  of  about  300  acres  adjacent  to  Terminal 

Island  with  a navigation  channel  and  mooring  slip  depth  of  ->4^  ft  mllw. 
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For  the  supertanker  facilities,  alternatives  may  range  from  dredging 
a 65-ft  channel  from  the  harbor  entrance  to  existing  tanker  facilities 
to  construction  of  three  new  supertanker  berths  and  dredging  of  a deeper 
navigation  channel.  New  landfill  would  be  added  in  Los  Angeles  Harbor 
with  material  dredged  from  the  new  channel. 

18.  In  Long  Beach  Harbor,  a proposed  development  by  the  Port  in- 
cludes deepening  of  the  main  navigation  channel  to  either  -65  or  -82  ft 
mllw,  addition  of  approximately  120  to  130  acres  of  landfill  southward 
of  Pier  J for  new  supertanker  and  container  ship  berths,  and  additional 
landfill  in  the  outer  harbor  between  the  Middle  breakwater  and  the  Navy 
Mole. 
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19-  An  example  of  an  early  proposed  development  plan  which  in- 
cluded a marine  development  by  the  City  of  Long  Beach  east  of  the  South- 
east Basin  is  shown  in  Figure  7.  The  development  plan  is  preliminary  but 
provides  an  indication  of  planned  harbor  improvements.  In  Long  Beach 
Harbor  several  revisions  including  reduction  in  the  size  of  the  landfill 
and  addition  of  a breakwater  have  been  proposed  adjacent  to  Pier  J.  Fu- 
ture development  of  the  harbor  will  undoubtedly  be  constructed  in  phases. 


Model  Description 


\ 


20.  The  Los  Angeles-Long  Beach  Harbors  hydraulic  model' (Figure  2) 
was  molded  in  cement  mortar  and  reproduced  to  scale  San  Pedro  Bay  and  a 
portion  of  the  Pacific  Ocean  surrounding  the  harbor  area.  The  model 
shoreline  extended  from  approximately  2 miles  northwest  of  Point  Fermin 
to  Huntington  Beach.  Underwater  contours  were  reproduced  out  to  the 


Figure  T.  Example  of  a proposed  development  plan  for  the 
Los  Angeles-Long  Beach  harbors 


( 
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-300  ft  rallw  contour,  and  sufficient  additional  offshore  area  was  in- 
cluded to  provide  room  for  wave  generators  and  the  tide  generator.  A 
general  view  of  the  model  is  shown  in  Figure  8.  Operation  of  the  tide 
generator  system  is  discussed  in  Reference  2. 

21.  The  model  scales  were  1:100  vertically  and  1:U00  horizon- 
tally. The  model  area  was  approximately  UU.OOO  sq  ft,  representing 
about  253  square  miles  of  prototype  area.  Depth  data  for  model  con- 
tours were  obtained  from  the  U.  S.  Coast  and  Geodetic  Survey  (now 
National  Ocean  Survey)  Charts  5101  and  51^7,  and  from  harbor  soundings 
provided  by  the  Ports.  Major  piers  and  wharves  were  reproduced  in  the 
model  with  l/l6-  and  1/32-in. -diam  brass  rods  used  to  simulate  pier 
piling.  The  bays  east  of  the  Harbor  such  as  Alamitos  Bay  and  Anaheim 
Bay  were  correctly  reproduced  in  plan  but  depths  were  averaged  in  the 
model  in  order  to  expedite  construction  and,  at  the  same  time,  permit 
proper  reproduction  of  approximate  tidal  volumes.  If  future  studies 
in  these  areas  are  needed,  the  actual  topography  can  be  installed  in 
the  model  with  relative  ease. 
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gure  8.  General  view  of  the  Los  Angeles 


PART  II:  MODEL  DESIGN  ANALYSIS 


Design  of  Model 

22.  For  analysis  of  wave-induced  harbor  oscillations  in  the 
Los  Angeles-Long  Beach  Harbors,  the  parameters  considered  during  model 
design  were  wave  refraction  and  diffraction,  viscous  friction  effects, 
wave  reflection,  wave  transmission,  and  wave  filters  and  wave  absorbers. 
Also  considered  were  the  absolute  dimensions  of  model  wave  heights  and 
periods  and  cost  of  model  construction. 

Wave  Refraction,  Diffraction,  and  Harbor  Response 

V 
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23.  Wave  refraction  is  the  process  by  which  the  direction  of 
waves  approaching  a coastal  area  at  an  angle  to  the  contours  is  changed. 
The  part  of  the  wave  advancing  in  shallower  water  moves  more  slowly 
than  that  part  advancing  in  deeper  water,  causing  the  wave  crest  to 
bend  toward  alignment  with  the  underwater  contours.  Dependent  upon  the 
shape  of  the  underwater  contours,  bending  of  the  wave  front  may  increase 
or  decrease  the  wave  height.  The  protection  provided  San  Pedro  Bay  by 
the  offshore  islands  and  the  broad  sill  west  of  San  Pedro  Bay  connect- 
ing San  Nicolas  and  Santa  Rosa  Islands  varies  with  wave  direction  and 
period.  However,  previous  refraction  studies^  have  shown  that  the  bay 
is  relatively  well  protected  from  long-period  wave  attack  except  for 
directions  between  south-southeast  and  south-southwest.  Wave  refraction 
analysis^  by  SEA  indicates  that  short-period  waves  (near  15  sec)  from 
the  south  are  focused  along  the  breakwater  near  the  Long  Beach  Harbor 
navigation  entrance  due  to  the  convex  shape  of  the  underwater  contours. 
SEA  wave  refractipn  diagrams^  from  the  south  for  longer  period  waves 
(30  sec  or  larger)  indicate  that  these  waves  can  propagate  into  the 
harbor  area  through  the  Gulf  of  Santa  Catalina  without  encountering  any 
effective  natural  barriers.  Wave  heights  seaward  of  the  harbor  area 
also  would  be  relatively  unaffected.  The  refraction  analysis  by  SEA 
also  included  other  directions  for  15-,  30- , 60-,  and  120-sec  waves,  and 
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the  results^  support  the  conclusion  that  San  Pedro  Bay  is  relatively 
well  protected  except  from  the  south- southeast  through  the  south- 
southwest  . 

13 

2**.  Diffraction  is  the  phenomenon  hy  which  energy  is  trans- 
mitted laterally  along  a wave  crest.  Near  an  island  or  breakwater,  the 
waves  propagate  into  the  sheltered  area  in  the  lee  of  the  structure. 

In  a distorted-scale  model,  refraction  and  diffraction  cannot  both  be 
properly  scaled  simultaneously.  Refractive  effects  depend  on  wave 
celerity  and  are  a function  of  h/L  (the  ratio  of  water  depth  h to  wave- 
length L).  Diffraction  is  a function  of  x/L  or  y/L  (the  ratio  of  hori- 
zontal distance  to  wavelength).  If  wavelength  is  scaled  according  to 
the  vertical  length  scale  in  a distorted-scale  model,  the  ratio  h/L 
will  be  in  exact  similitude  and  consequently  refraction  similitude  will 
be  achieved.  If  wavelength  is  scaled  according  to  the  horizontal 
length  scale  in  a distorted-scale  model,  the  ratio  x/L  or  y/L  will  be 
in  exact  similitude  and  diffraction  similitude  will  be  achieved.  In 
the  Los  Angeles-Long  Beach  model  study,  the  most  important  design  con- 
straint is  that  similitude  of  mode  shapes  and  resonant  frequencies  of 
oscillation  be  precisely  achieved.  The  governing  Helmholtz  equation^ 
for  harbor  oscillations  is 


where 

x,y,z  = axes  of  a rectangular  coordinate  system  fixed  at  the 
mean  water  surface 

n = local  surface  elevation 

o = angular  frequency 

g = acceleration  due  to  gravity 

Since  the  same  equation  applies  in  model  and  prototype.  Equation  1 may 
be  written  as 
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+ 


+ n 

r r 


np  = 0 


where  the  subscript  p represents  the  prototype  and  r , the  scale 
ratio  of  model  to  prototype.  From  inspectional  analysis,  the  coeffi- 
cients of  Equation  2 must  be  equal,  or 


after  dividing  through  by  . Equation  3 indicates  that  a hydraulic 
model  may  be  distorted  for  proper  simulation  of  harbor  resonant  oscilla- 
tion frequencies.  The  angular  frequency  may  be  written  in  terms  of 
wavelength  and  water  depth  and  Equation  2 indicates  that  the  wavelength 
must  be  scaled  by  the  horizontal  scale  ratio.  Wave  height,  however,  may 
be  scaled  by  the  vertical  scale  ratio. 

25.  Therefore,  using  the  horizontal  scale  ratio  L for  wave- 

r 

length,  diffraction  and  harbor  oscillation  frequencies  will  be  in 
similitude  but  refraction  will  have  a scale  effect  due  to  model  dis- 
tortion. For  shallow-water  waves  where  the  wave  celerity  is  governed 
by  the  local  depth,  model  distortion  will  have  little  effect  on  re- 
fraction. The  equation  for  wave  celerity  C from  small  amplitude  wave 
theory  is 

. .1/2 

C = tanh  khj  (1) 

where  k is  the  wave  number,  2ir/L  . As  the  wave  period  increases, 
tanh  kh  approaches  kh  and,  for  shallow- water  waves,  the  celepity 
becomes 


C = 
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Consequently,  for  shallov-water  waves,  model  distortion  has  a negligible 
effect  on  wave  refraction.  The  ratio  of  celerity  calculated  from  Equa- 
tion ^ to  the  shallow-water  approximation  of  Equation  5 is  shown  in 
Eigure  9 as  a function  of  wave  period  for  water  depths  of  39,  50,  100, 
and  300  ft  (for  both  the  prototype  and  the  distorted-scale  model).  As 
indicated  by  Eigure  9,  the  distortion  scale  effect  on  wave  celerity  is 
negligible  (<1.3^  at  300  ft)  for  T ^ 120  sec.  At  depths  less  than 
300  ft,  the  wave  period  above  which  the  wave  celerity  remains  relatively 
constant  (i.e.  C -+  /^)  decreases.  Tne  shape  of  the  offshore  topog- 
raphy also  influences  wave  refraction  patterns,  and  if  the  wave  front 
propagates  approximately  normal  to  the  offshore  contours,  refraction 
patterns  may  be  similar  at  periods  much  less  than  the  lower  limit  in- 
dicated for  1.3  percent  accuracy  by  the  shallow-water  celerity 
approximation. 

26.  Due  to  the  model  distortion,  the  time  scale  ratio  T is 

r 

affected  by  the  choice  of  vertical  or  horizontal  scale  for  wavelength. 
The  angular  wave  frequency  is 


o = ^ = (g^  tanh  kh)^'^^ 


(6) 


Using  the  horizontal  scale  for  wavelength,  the  time  scale  ratio  can 
then  be  written  as 


(7) 


where  the  subscript  m refers  to  model  quantities.  The  equation  can 
also  be  written  as 


T 

r 


(Ta) 
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Figure  9*  Ratio  of  wave  celerity  calculated  from  the  shallow-water  approximation 

for  depths  of  39,  50,  100,  and  300  ft 


where  n is  the  model  distortion.  As  tanh  kh  approaches  kh  , the 
time  scale  ratio  may  be  approximated  by 


T 

r 


(7b) 


For  an  average  depth  of  39  ft,  the  prototype  wave  period  calculated 
from  Equation  7b  is  within  1 percent  of  the  period  calculated  from 
Equation  7 for  T ^ 85  sec  . In  a distorted-scale  model  where  refrac- 
tion is  the  most  significant  effect  and  wavelength  is  scaled  by  the 
vertical  scale,  the  time  scale  ratio  is: 


T = 
r 


(7c) 


27.  In  summary,  for  analysis  of  wave-induced  harbor  oscillations, 
the  preceding  discussion  of  simulating  refraction,  diffraction,  and 
resonant  frequencies  of  oscillation  in  a distorted  model  indicates  that: 

a.  Wavelengths  in  the  model  should  be  scaled  by  the  hori- 
zontal scale  ratio. 

b.  The  time  scale  ratio  from  model  to  prototype  is  depen- 
dent on  wave  period  and  water  depth  and  is  given  by 
Equation  7. 

c_.  The  harbor  is  relatively  well  protected  from  long-period 
wave  attack  except  from  the  south-southeast  through 
south-southwest . 

Refraction  in  the  distorted-scale  model  and  prototype 
will  be  in  similitude  (to  a very  close  approximation) 
for  longer  wave  periods  (say,  greater  than  75  sec). 
However,  at  shorter  wave  periods,  there  may  be  some 
refraction  scale  effects  that  require  adjustment  in 
either  the  wave  generator  position  and/or  wave  height 
in  order  to  ensure  similitude  of  wave  refraction. 


Results  of  a detailed  wave  refraction  analysis  are  given  in  the  next 
sections  for  distorted-  and  undlstorted-scale  conditions. 

Wave  Refraction  Analysis 

28.  A detailed  wave  refraction  analysis  was  conducted  for  the 
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offshore  area  using  an  analysis  procedure  that  calculates  wave  orthog- 
onals  with  the  usual  assumption  that  wave  energy  is  not  transmitted 
along  the  wave  front  perpendicular  to  the  wave  orthogonal . Wave  height 
H is  given  by 


H = K K H. 
r s 1 


(8) 


where 

K = refraction  coefficient 
r 

= shoaling  coefficient 
= initial  wave  height 

Bottom  reflection  is  neglected  since  this  effect  is  normally  not  sig- 
nificant for  small  slopes.  Bottom  friction  and  percolation  also  are 
neglected^^  since  their  effect  is  minor  in  the  prototype  for  small 
amplitude,  long-period  waves.  As  a part  of  the  procedure,  a second- 
degree  polynomial  is  fitted  by  a least-squares  technique  to  12  points 
in  a rectangular  grid  with  a square  mesh  surrounding  the  area  through 
which  the  orthogonal  is  propagating. 

29.  The  accuracy  of  the  polynomial  representation  of  complex 
offshore  topography  increases  as  the  number  of  points  in  the  rectangular 
grid  increases.  Two  rectangular  grids  were  used  for  prototype  condi- 
tions in  the  wave  refraction  analysis:  Grid  1 covers  a large  area  with 

a wide  grid  spacing;  Grid  2 covers  a small  area  adjacent  to  the  harbor 
with  a more  detailed  grid;  and  the  boundaries  of  the  two  grids  are  shown 
in  Figure  10.  The  wave  refraction  data  from  Grid  1 are  used  for  the 
initial  wave  orthogonal  data  in  Grid  2.  Grid  point  spacing  and  grid 
size  were: 

Grid  No.  Point  Spacing,  ft  Grid  Size,  ft 

1 9TU0  360,380  X 1*18,820 

2 1675  110,550  X 187,600 

Depth  data  at  each  grid  point  were  obtained  from  the  U.  S.  Coast  and 
Geodetic  Survey  (now  National  Ocean  Survey)  Chart  5101  for  Grid  1 and 
Chart  51^2  for  Grid  2.  Grid  2 data  also  were  used  for  the  physical 
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Figure  10-  Location  of  wave  refraction  Grids  1 and  2 


model  wave  refraction  analysis  with  depths  limited  to  a maximum  of 
300  ft  mllw. 

30.  Wave  refraction  diagrams  were  calculated  for  15,  20,  30,  ^*5, 
60,  75,  90,  120,  l80,  2l<0,  300,  and  360  sec  using  Grid  1 for  waves  ap- 
proaching from  the  south.  Plates  1-1<  show  the  refraction  diagrams  for 
wave  periods  of  15,  30,  60,  and  300  sec.  Refraction  patterns  for  the 
remaining  periods  are  similar.  The  diagrams  show  the  effect  of  Santa 
Catalina  Island  on  wave  refraction  at  the  longer  wave  periods.  Strong 
convergence  of  wave  energy  seaward  of  the  Queens  Gate  breakwater  en- 
trance is  vividly  illustrated.  The  model  limits  also  are  shown  in 
Plates  1-^*  and  have  been  selected  to  ensure  that  the  convergent  zone 
is  included  within  the  model  area.  Within  the  strong  convergent  zone, 
diffraction  effects  not  accounted  for  by  normal  wave  refraction  analyses 
become  quite  significant  ’ and  refraction  coefficients  are  no  longer 
valid.  The  convergent  zone,  although  not  as  pronounced,  exists  for  a 
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wave  period  of  15  sec  as  well  as  longer  periods  and  is  a result  of  the 
convex  bottom  topography  and  submarine  canyons.  The  -300  ft  mllw  con- 
tour encloses  the  convergence  zone. 

31.  After  selecting  the  model  limits  based  on  the  preceding  wave 
refraction  data,  a more  detailed  analysis  of  wave  refraction  was  neces- 
sary within  the  model  limits.  Using  Grid  2 (Figure  10),  new  wave  orthog 
onals  were  started  (using  initial  conditions  based  on  the  data  from 
Grid  l)  for  selected  wave  fronts.  Wave  refraction  diagrams  again  were 
calculated  for  15- , 20- , 30- , U5-,  60-,  75-,  90- , 120- , l80-,  2^0-, 

300- , and  360-sec  wave  periods.  Examples  of  these  diagrams  for  30,  60, 
l80,  and  360  sec  are  shown  in  Plates  5-8.  The  orthogonal  convergent 
zone  is  shown  in  more  detail  in  these  diagrams.  At  the  lower  periods 
(approximately  15  to  30  sec)  the  refraction  diagrams  indicate  that  the 
convergent  zone  is  approximately  focused  on  Queens  Gate.  As  the  period 
increases,  however,  the  focus  of  the  convergent  zone  shifts  slightly 
eastward  along  the  Long  Beach  breakwater.  Provided  wave  periods  are 
greater  than  approximately  75  sec,  refraction  patterns  are  similar  near 
the  convergent  zone.  The  wave  orthogonals  do  not  directly  overlay  since 
the  orthogonals  were  started  from  varying  locations  along  the  wave  front 
selected  for  each  wave  period  due  to  refraction  occurring  seaward  of  the 
initial  wave  front  used  for  the  detailed  analysis.  Similarity  of  re- 
fraction patterns  (for  T ^ 75  sec)  is  primarily  due  to  the  effect  of 
depth  on  wave  celerity  discussed  in  paragraph  25. 

32.  The  shape  of  the  wave  front  for  15-,  20-,  30-,  I+5-,  60-,  75-, 
and  360-sec  wave  periods  immediately  after  the  wave  front  enters  within 
the  model  limits  and  approaches  the  -300  ft  mllw  contour  is  shown  in 
Plates  9-11.  The  15-  and  20-sec  wave  fronts  are  relatively  straight  but 
the  30-  and  1+5-sec  wave  fronts  are  starting  to  curve  as  the  refractive 
effects  increase.  For  longer  wave  periods,  shape  of  the  wave  fronts  is 
similar  near  the  -300  ft  mllw  contour.  However,  near  the  east  and  west 
model  limits  wave  front  shape  varies  with  wave  period.  This  variation 
is  due  to  the  effect  of  Santa  Catalina  Island,  the  two  submerged  canyons 
and  nearshore  topography  at  the  east  and  west  model  limits.  The  varia- 
tion also  can  be  seen  in  the  large  grid  wave  refraction  diagrams  shown 
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in  Plates  1-U.  In  general,  bathymetry  near  Santa  Catalina  Island  causes 
the  wave  orthogonals  passing  the  island  to  bend  northwestward  as  the 
period  increases.  Nearshore,  the  orthogonals  tend  to  bend  perpendicular 
to  the  offshore  topography  as  the  period  increases.  The  complexity  of 
these  effects  precludes  a consistent  trend  in  wave  front  shape  (as  a 
function  of  period)  near  the  model  limits.  It  will  be  illustrated  later 
that  differences  in  wave  front  location  at  longer  wave  periods  (in  the 
prototype)  seem  to  be  balanced  by  the  effects  of  model  distortion  and 
limitation  of  model  depths  to  -300  ft  mllw,  and  consequently,  a single 
wave  front  position  can  be  used  in  the  model  for  longer  periods. 

33.  The  effect  of  model  scale  distortion  on  wave  refraction  was 
initially  investigated  using  Grid  2 with  the  grid  point  distance  scaled 
by  the  horizontal  scale  and  depths  by  the  vertical  scale  for  the  follow- 


ing  conditions: 
Condition  No. 

Distortion  Ratio 

Vertical  Scale 

Horizontal  Scale 

1 

1:2 

1:100 

1:200 

2 

1:3 

1 ; 100 

1 : 300 

3 

1:1* 

1:100 

1:1*00 

As  indicated  above,  model  distortion  ratios  of  1:2,  1:3,  and  1:1*  were 
considered.  Change  in  model  scale  with  the  same  distortion  ratio  pro- 
duces identical  refraction  patterns  and  is  not  considered.  The  1:100 
vertical  scale  was  near  the  minimum  scale  ratio  necessary  for  measure- 
ment of  small  amplitude  waves.  Typical  examples  of  distorted-scale 
wave  refraction  diagrams  are  shown  in  Plates  12-23  for  conditions  1 to 
3 at  30- , 60-,  l80-,  and  360-sec  wave  periods.  The  distorted-scale 
refraction  analysis  also  included  the  15-,  20-,  1*5-,  75-,  90-,  120-, 
2l*0-,  and  300-sec  periods.  The  initial  data  for  orthogonal  azimuth  and 
wave  front  starting  location  were  taken  from  the  results  of  the 
undistorted-scale  (prototype)  Grid  2 wave  refraction  analysis.  Com- 
parison of  the  distorted-scale  wave  refraction  diagrams  indicates  that 
wave  refraction  patterns  for  the  1:3  and  1:1*  distortion  ratios  are 
similar  for  all  wave  periods.  However,  at  15  and  30  sec,  the  wave 
refraction  pattern  for  the  1:2  distortion  ratio  varies  significantly 
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from  the  two  larger  distortion  ratios.  At  and  above  60  sec,  the  refrac- 
tion patterns  are  similar  for  all  three  distortion  ratios.  Also,  com- 
parison of  the  refraction  diagrams  for  the  prototype  with  the  corre- 
sponding refraction  diagrams  for  distorted-scale  conditions  at  15- , 

30- , 1*5-,  and  60  -sec  wave  periods  indicates  that  the  prototype  and 
distorted-scale  refraction  will  be  similar  above  approximately  60  sec; 
but  near  60  sec  and  below,  the  location  of  the  wave  front  for  distorted- 
scale  conditions  is  affected  by  the  distortion  and  must  be  adjusted  in 
order  for  model  wave  refraction  to  simulate  the  prototype  refraction. 

3I*.  Due  to  the  results  of  the  preceding  analysis  of  distortion 
effects  on  wave  refi’action,  a model  distortion  ratio  of  1:1*  vas  selected. 
The  1:1*  distortion  ratio  minimizes  the  model  area  and  does  not  signifi- 
cantly increase  the  period  range  over  which  adjustment  of  the  initial 
model  wave  front  is  necessary  for  refraction  similitude.  It  was  clearly 
illustrated  that  there  was  a negligible  difference  in  refraction  scale 
effects  when  changing  the  distortion  factor  from  3 to  1*;  however,  there 
was  a difference  for  a distortion  ratio  of  2.  Thus,  the  only  question 
is  whether  to  choose  a distortion  of  2 or  1*.  The  cost  reduction 
achieved  by  the  large  reduction  in  model  area  (a  factor  of  1*)  using  a 
1*  distortion  ratio  greatly  outweighs  the  relatively  modest  technical 
advantage  of  slightly  smaller  refraction  scale  effects.  This  is  es- 
pecially true  when  the  refraction  scale  effect  can  be  corrected  (as  is 
the  case  here)  by  a change  in  the  model  wave  front  position.  Since  the 
wave  refraction  convergent  zone  is  inclosed  within  the  -300  ft  mllw 
contour,  the  underwater  contours  were  reproduced  out  to  this  contour. 

The  submerged  canyons  near  the  east  and  west  model  limits  are  not  re- 
produced (below  the  -300  mllw  contour)  and  adjustment  of  the  initial 
model  wave  front  to  compensate  for  the  decreased  depths  in  the  canyon 
areas  is  necessary. 

35"  Adjustment  of  the  initial  model  wave  front  position  was  ac- 
complished by  an  iterative  procedure.  A wave  front  location  and  wave 
orthogonal  azimuth  were  selected,  a refraction  diagram  calculated,  a 
comparison  made  with  the  refraction  diagram  calculated  using  Grid  2 for 
the  same  period,  and  the  process  was  repeated  until  approximate 
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similitude  of  wave  fronts  and  orthogonal  azimuths  was  achieved  seaward 
of  the  breakwater  and  adjacent  to  the  convergent  zone.  The  initial 
prototype  and  adjusted  wave  front  locations  are  shown  in  Plates  9-11 
for  15-,  20-,  30-,  60-,  75-,  and  360  -sec  wave  periods.  Also  shown 

in  Plates  9-H  are  the  comparisons  of  wave  front  positions  for  each  of 
the  wave  periods  seaward  of  the  harbor  breakwater  after  the  convergent 
zone  has  formed. 

36.  In  general,  the  prototype  and  model  wave  fronts  near  the 
breakwater  compare  closely  at  the  lower  wave  periods  and  reasonably  well 
as  the  period  increases.  For  the  longer  periods  shown  in  Plates  10  and 
11 , the  model  wave  fronts  tend  .to  diverge  slightly  from  the  prototype 
near  either  end.  While  it  may  appear  that  minor  shifting  of  the  initial 
model  wave  front  would  result  in  more  exact  agreement,  the  prototype 
wave  orthogonal  azimuths  along  the  wave  front  also  must  be  simulated 
and  the  intial  wave  front  shown  represents  a compromise  between  the 
requirements  that  both  wave  front  position  and  orthogonal  azimuth  should 
be  in  similitude.  Examples  of  the  agreement  between  prototype  and  model 
orthogonals  are  given  in  Tables  1 and  2 for  20-  and  T5-sec  wave  periods. 
As  indicated  by  the  azimuth  data,  the  maximum  variation  of  the  model  az- 
imuth from  the  prototype  for  the  20-sec  wave  period  is  9.1  degrees  with 
an  average  variation  of  3.9  degrees.  Four  of  the  ten  20-sec  orthogonals 
have  a model  azimuth  less  than  the  prototype.  For  the  75-sec  period, 
the  maximum  azimuth  variation  is  8.1  degrees.  The  variation  for  the 
remaining  l.T  Oithogonals  is  less  than  5 degrees  and  the  average  varia- 
tion is  2.8  degrees.  Five  of  the  thirteen  75-sec  orthogonals  have  a 
model  azimuth  less  than  the  prototype. 

37.  The  adjusted  wave  front  location  shown  in  Plates  9-11  is 
composed  of  straight  line  segments  corresponding  to  the  length  of  each 
unit  of  the  wave  generator.  The  straight-line  segments  also  affect  the 
accuracy  within  which  the  wave  front  location  and  orthogonal  azimuth 
can  be  matched. 

38.  As  indicated  in  Plates  9-11,  the  initial  model  wave  front 
locations  for  15-,  20-,  and  30-sec  wave  periods  show  little  variation, 
especially  in  the  central  sections  of  the  wave  front  near  the  -300  ft 
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contour;  thus,  in  the  15-  to  30-sec  range,  an  average  of  the  three  ini- 
tial wave  front  locations  was  selected  for  the  initial  model  wave  front 
location.  Similarly,  an  average  wave  front  was  selected  for  the  ^45-  to 
60-sec  range.  These  average  initial  wave  fronts  are  shown  in  Plate  2h 
along  with  the  initial  model  wave  front  for  the  T5-sec  and  above  wave 
period  range.  As  shown  in  Plate  2h , the  average  of  the  initial  model 
wave  front  for  the  I45-  to  60-sec  range  is  shifted  significantly  from 
the  15-  to  30-sec  average  initial  wave  ftront  along  the  western  limit  of 
the  wave  front.  Additional  model  wave  refraction  analysis  was  necessary 
to  determine  the  initial  model  wave  front  to  be  used  in  the  transition 
from  30  to  i*5  sec.  The  wave  front  selected  as  a result  of  the  additional 
analysis  for  35  sec  also  is  shown  in  Plate  2h  and  provides  a ..mooth 
transition  from  the  15-  to  30-sec  range  to  the  !+5-  to  60-sec  range. 

The  additional  analysis  indicated  the  average  ^5-  to  60-sec  model  wave 
front  was  acceptable  for  a UO-sec  period. 

39-  Wave  heights  at  each  of  the  wave  generators  also  were  calcu- 
lated in  the  refraction  analysis.  The  prototype  wave  height  was  cal- 
culated using  Equation  8 at  a location  seaward  of  the  convergent  zone 
and  within  the  -300  ft  contour.  Then,  normalized  wave  heights  along  the 
wave  front  at  the  wave  generator  were  calculated  using  model  refraction 
coefficient  data.  The  normalized  wave  heights  for  each  of  the  1^4  wave 
generator  units  at  15,  20,  30,  ^45,  60,  75,  and  90  to  360  sec  are  shown 
in  Table  3.  The  shoaling  coefficient  has  not  been  used  in  developing 
Table  3.  There  is  a shoaling  coefficient  scale  effect  in  the  model 
(due  to  the  scale  distortion  and  method  of  scaling)  which  will  be  ac- 
counted for  in  analyzing  the  model  data;  however,  shoaling  does  not 
affect  the  normalized  wave  height  at  the  wave  generator. 

^0.  Normalized  wave  heights  tend  to  reach  a maximum  near  wave 
generator  units  10,  11,  and  12  in  the  15-  to  30-sec  period  range;  but 
at  U5  sec  and  above,  the  maximum  normalized  value  occurs  near  wave 
generator  unit  12.  Location  of  wave  generator  units  in  the  15-  to 
30-sec  period  range  are  shown  in  Figure  11.  The  shift  in  the  maximum 
normalized  wave  height  to  wave  generator  unit  12  is  probably  due  to  the 
increased  northwestward  refraction  near  Santa  Catalina  Island  discussed 
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Figure  11.  Location  of  wave  generator  units 
1 to  lU  in  15-  to  30-sec  period  range 

in  paragraph  31.  The  minimum  normalized  wave  height  along  the  initial 
model  wave  front  occurs  at  wave  generator  unit  2 and  is  probably  due  to 
the  effect  of  the  submerged  canyon  on  the  western  side  of  the  harbor 
area  on  wave  refraction.  The  normalized  wave  height  data  also  indicate 
that  the  height  of  the  wave  front  propagating  into  the  Los  Angeles 
Harbor  area  is  less  than  the  height  of  the  wave  front  propagating  into 
the  Long  Beach  Harbor  area.  For  a 15-sec  wave  period,  the  refraction 
analysis  indicates  that  for  an  initial  wave  front,  the  wave  height  at 
Angels  Gate  should  be  approximately  Uo  percent  of  the  wave  height  at 
Queens  Gate.  If  the  variation  in  wave  height  along  the  wave  front  is 
not  simulated,  the  resonant  frequencies  in  the  harbor  area  will  not  be 
affected  but  resonant  response  Eimplitudes  and  the  wave  energy  distribu- 
tion could  be  affected.  Thus,  it  was  planned  to  simulate  the  variable 
height  along  the  wave  front. 
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Ul.  A brief  summary  of  the  conclusions  from  the  wave  refraction 
analysis  are: 

Santa  Catalina  Island  has  a significant  effect  on  wave 
refraction,  particularly  at  the  longer  wave  periods. 

b.  A convergent  zone  forms  seaward  of  the  Middle  and  San 
Pedro  breakwaters  over  the  15-  to  360-sec  period  range 
included  in  the  analysis. 

c_.  Refraction  patterns  were  similar  at  and  above  a wave 
period  of  about  75  sec  near  the  convergent  zone. 

For  scale  distortion  ratios  of  1:3  and  1:^,  the  calcu- 
lated refraction  patterns  were  similar  over  the  15-  to 
360-sec  period  range,  but  for  the  1:2  distortion  ratio, 
the  calculated  refraction  was  significantly  different  at 
lower  periods. 

£.  Near  a wave  period  of  60  sec  and  below,  the  calculated 
refraction  patterns  for  distorted  scale  conditions 
changed  significantly  from  the  calculated  prototype  re- 
fraction patterns,  and  adjustment  of  the  initial  wave 
front  in  the  model  was  necessary. 

f^.  A model  distortion  ratio  of  l:l4  and  a vertical  scale 

ratio  of  1:100  were  selected  to  minimize  model  area  and 
to  provide  a vertical  scale  ratio  where  accurate  model 
measurements  could  be  assured. 

The  position  of  the  initial  model  wave  front  must  be  ad- 
justed to  compensate  for  distortion  scale  effects  and 
for  refraction  caused  by  bathymetry  at  depths  greater 
than  -300  ft  mllw  (limit  of  model  bathymetry). 

Normalized  wave-height  variation  along  the  prototype  wave 
front  is  significant  and  should  be  reproduced  in  the 
model  wave  tests. 

Model  Wave-Height  Attenuation 

h2.  Waves  in  a model  may  be  attenuated  by  surface  tension,  inter- 
nal friction,  and  friction  in  the  bottom  boundary  layer.  In  the  wave 
period  range  to  be  simulated  in  this  study,  surface  tension  effects  in 
crease  as  the  wave  period  decreases;  and  for  a model  wave  period  of 

0.5  sec,  near  the  minimum  to  be  generated  in  the  model,  the  effect  of 
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surface  tension  on  wave  celerity  is  less  than  0.1  percent.  Internal 
friction  effects  decrease  as  the  wave  period  increases,  and  for  the 
0.5-sec  model  wave  period,  the  wave-height  attenuation  from  wave 
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generator  unit  9 to  the  Long  Beach  breakwater  is  less  than  3.2  percent 
For  a 1-sec  wave  period,  the  height  attenuation  is  less  than  0.2  per- 
cent. Using  an  average  harbor  depth  of  39  ft,  the  0.9-  and  1-sec  wave 
periods  represent  prototype  periods  of  lU.L  sec  and  37.0  sec,  respec- 
tively. The  wave-height  attenuation  due  to  internal  friction  was  cal- 
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culated  from 


^t  -8TT^Vt/L^ 

17=  ^ 


where  H is  the  wave  height  at  time  t and  v is  the  kinematic  vis- 

^ _c 

cosity.  The  kinematic  viscosity  was  assumed  to  be  1.059  (lO  for  T0°F 

i^3.  The  attenuation  of  wave  height  for  a progressive  oscillatory 

wave  with  distance  x due  to  friction  in  the  bottom  boundary  layer  may 
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be  calculated  from  Keulegan's  equation 


where 


H, 


-ax 

e 


( ■ V,  . ^4nh  \ 

(sinh  — + —} 


(10) 


(11) 


and  T is  the  wave  period.  In  developing  Equation  11,  the  width  of  the 
wave  crest  is  assumed  to  be  many  times  greater  than  the  depth  in  a rect- 
angular channel  of  uniform  cross  section.  The  attenuation  factor  at 
various  periods  for  two  locations  is  tabulated  below.  The  first  attenu- 
ation factor  column  is  for  a wave  traveling  from  wave  generator  unit  1 
in  the  ^5-  to  60-sec  position  to  the  San  Pedro  breakwater  which  has  a 
minimum  wave  travel  distance.  The  second  attenuation  factor  column  is 
for  a wave  traveling  from  wave  generator  unit  9 to  the  Long  Beach  break- 
water and  is  representative  of  the  maximum  wave  travel  distance  before 
entering  the  harbor. 
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Attenuation  Factor 


H 


i 


San  Pedro 

Long  Beach 

Period,  sec 

Breakwater 

Breakwater 

15 

0.99 

0.99 

20 

0.99 

0.97 

30 

0.98 

0.9*+ 

1*5 

0.97 

0.92 

60 

0.97 

0.93 

75 

0.97 

0.93 

90 

0.98 

0.9'4 

120 

0.98 

0.95 

180 

0.99 

0.97 

2l*0 

0.99 

0.98 

300 

0.99 

0.98 

360 

0.99 

0.99 

1*1*.  The  maximum  variation  in  the  attenuation  factor  for  the 

two  travel  distances  is  0.05, 

indicating  that  variation  in  travel  dis- 

tance  for  the  present  case  has 

a relatively 

small  effect  on  wave-height 

attenuation.  The  wave  height 

may  be  varied 

along  the  generated  wave 

front  to  correct  for  the  small 

variation  in 

wave-height  attenuation  or. 

when  comparing  model  wave  heights  with  prototype  heights  directly,  the 


observed  model  wave  heights  can  be  increased  to  reflect  attenuation  in 
the  model. 


Wave  Reflections 


1*5.  Wave  reflections  may  occur  from  the  underwater  topography 
seaward  of  the  harbor,  from  the  breakwaters,  and  from  the  harbor  side 
slopes.  In  the  prototype,  wave  reflection  from  the  underwater  topog- 
raphy and  the  breakwater  would  propagate  seaward;  but  in  the  model, 
waves  reflected  from  these  sources  may  also  re-reflect  from  the  model 

wave  generator  and  possibly  affect  the  incident  waves.  Reflections  from 

21 

the  model  side  slopes  within  the  harbor  may  affect  the  amplitude  of 
harbor  oscillations  and  should  be  correctly  simulated.  The  wave  reflec- 
tion coefficient  is  the  ratio  of  the  reflected  wave  height  to 

the  incident  wave  height  or 
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(12) 


H 

1 

The  reflection  coefficient  varies  from  1.0  for  total  reflection  to  0.0 
for  no  reflection.  In  general,  the  magnitude  of  the  reflection  coef- 
ficient will  depend  on  geometry,  composition  of  a structure,  and  inci- 
dent wave  characteristics  such  as  wave  steepness  n/L  , (n  is  the  local 
water-surface  elevation)  and  relative  depth  h/L  near  the  structure. 

46.  Reflection  from  the  offshore  contours  was  estimated  using  a 
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theory  developed  hy  Rosseau  which  provides  an  analytical  solution  for 
the  reflection  coefficient.  The  development  of  the  theory  also  is  de- 
scribed in  Reference  l6.  Rosseau' s theory  requires  no  assumptions  con- 
cerning wave  form  other  than  linear  small  amplitude  progressive  waves 
propagating  from  +”  in  a canal  of  variable  depth  as  shown  in  Figure  12, 
and  that  the  bottom  profile  may  be  described  by  a specified  function 
such  that  the  depth  approaches  as  y goes  to  and  h^  as  y 

goes  to  +“  . The  solution  is  not  valid  for  a discontinuous  change  in 
slope;  however,  the  Rosseau  contour  can  approach  the  shape  of  a discon- 
tinuous sloped  step  in  bottom  contours.  An  incompressible  nonviscous 
fluid  is  assumed  and  the  solution  is  obtained  using  a conformal  mapping 
technique. 

47.  In  the  bottom  contour  reflection  analysis,  the  following 
three  cases  were  considered  for  a typical  model  depth  variation  with  a 
length  of  50  ft  and  a change  in  depth  of  1.5  ft: 
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Case  No. 

Depth  hj^  , ft 

Depth  h^  , ft 

1 

2.0 

0.5 

2 

2.5 

1.0 

3 

3.5 

2.0 

These  depths  are  representative  of  steep  model-depth  changes  seaward  of 
the  harbor  breakwater.  In  Figure  13,  a comparison  of  the  Rosseau  con- 
tour with  the  bottom  contour  for  Case  1 shows  that  the  Rosseau  contour 
closely  fits  the  bottom  slope.  The  comparisons  for  Cases  2 and  3 are 
similar . 
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Figure  13.  Rosseau  contour  for  Case  1 


U8.  Reflection  coefficients  for  various  model  wave  periods  calcu- 
lated for  the  three  cases  are  given  in  the  following  tabulation: 


Model  Wave 
Period,  sec 

Case  1 

Case  2 

Case  ' 

0.75 

<0.01 

<0.01 

<0.01 

5.00 

0.20 

0.15 

0.10 

10.00 

0.29 

0.  21 

0.13 

Using  an  average  depth  in  the  harbor  area  of  39  ft,  the  wave  periods  in 
the  table  correspond  to  26. 0-,  199.9-,  and  l^OO. 0-sec  prototype  wave  pe- 
riods. The  reflection  coefficients  range  from  essentially  0 in  the 
lower  wave  period  range  to  0.29  for  Case  1.  The  reflection  coefficients 
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for  all  three  cases,  however,  indicate  that  wave  filters  in  front  of  the 
wave  genersitor  will  be  required  at  longer  periods  to  attenuate  the 
height  of  the  reflected  wave  from  changes  in  the  offshore  topography. 

The  period  range  over  which  wave  filters  will  be  required  is  also  af- 
fected by  reflection  from  the  harbor  breakwaters  and  model  sidewalls. 
Many  sections  of  the  steep  model  depth  changes  are  oriented  such  that 
the  wave  is  reflected  toward  the  model's  perimeter  walls,  and  wave  ab- 
sorbers will  be  required  along  the  walls  to  attenuate  these  reflections. 

i+9.  Wave  reflections  from  the  distorted-scale  harbor  breakwater 
also  may  require  attenuation  by  the  wave  filters.  Two-dimensional  wave 
flume  tests  conducted  in  the  20-  to  100-sec  prototype  wave  period  range 
(described  in  paragraph  56)  to  determine  the  reflection  coefficient  from 
the  Middle  and  Long  Beach  model  breakwaters  indicate  the  reflection  co- 
efficient varies  from  approximately  0. U to  0.6.  Two-dimensional  fliime 
tests  were  not  conducted  above  100  sec  due  to  flume  length  limitations. 
The  model  wave  reflection  coefficients  determined  from  the  two- 
dimensional  tests  are  given  in  Table  U for  20-,  30-,  ^tO-,  50-,  60-,  T9-, 
80-,  90- , and  100-sec  prototype  wave  periods.  The  reflection  coeffi- 
cients given  in  Table  U are  for  prototype  wave  heights  of  approximately 
1 ft  and  for  approximately  3 to  U ft.  Design  of  the  model  breakwater  is 
discussed  in  more  detail  in  the  following  breakwater  transmission  sec- 
tion. Two-dimensional  tests  were  not  conducted  for  the  San  Pedro  break- 
water section;  however,  the  reflection  coefficient  should  decrease  with- 
out an  impervious  core. 

50.  To  evaluate  the  effect  of  side-slope  distortion  in  the  harbor 
area,  two-dimensional  tests  in  a 2-ft-wide  flume  with  a flapper-type 
wave  generator  were  conducted  for  a typical  prototype  side  slope  (Fig- 
ure ll<)  with  a vertical  wall  down  to  -5  ft  and  a l-on-1,5  slope  from 
-5  to  -50  ft  at  a 1:25  scale.  Still-water  level  during  the  tests  was  at 
+5  ft.  Wave  data  from  a moving  wave  gage  were  recorded  on  strip  charts 
during  the  tests.  Reflection  coefficients  were  calculated  from 
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where 


H = wave  height  at  the  antinode  of  the  wave  envelope  recorded  hy 
the  moving  wave  gage 

H = height  at  the  node 
n 


Figure  l^*.  Typical  prototype  and  model  side-slope  section 

51.  The  reflection  coefficients  calculated  using  Equation  13  for 
various  periods  are  as  follows. 


Prototype  Wave 
Period,  sec 

Prototype  Wave 
Height,  ft 

Reflection 

Coefficient 

20 

0.93 

O.8I4 

30 

0.96 

0,80 

UO 

1.05 

0,78 

50 

0.97 

0.78 

60 

I.2I+ 

0.75 

In  some  cases,  harmonics  observed  on  the  strip  chart  wave  record  made 

determination  of  the  wave  height  at  the  node  H difficult.  Generally, 
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in  tie  presence  of  harmonics  at  the  node,  the  reflection  coefficient 

was  actually  higher  than  that  calculated  from  Equation  13. 

52.  If  the  typical  side-slope  section  is  considered  as  a smooth 

slope  with  a ’-on-1.5  slope,  the  theoretical  reflection  coefficient,  cal- 
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culated  from  the  Miche  theory  for  wave  reflection  from  a plane  sloping 
beach,  is  1.0  for  the  periods  of  20  sec  or  longer  and  wave  steepness 
values  H/L  Included  in  the  experimental  test  series.  For  a flat  side 
slope,  the  reflection  coefficient  is  usually  modified  by  an  empirical 
coefficient  (<l)  which  is  constant  for  a given  structure.  For  a flat 


side  slope,  a value  of  0.9  can  be  used  and  the  resulting  theoretical 
reflection  coefficient  for  small  wave  steepness  (as  in  this  case)  is  0.9 
at  periods  of  20  sec  and  longer.  Allowing  for  experimental  error,  the 
calculated  reflection  data  from  the  two-dimensional  wave  tests  for  the 
typical  side  slope  and  the  theoretical  reflection  coefficient  for  a flat 
side  slope  are  in  good  agreement  and  indicate  that  the  typical  side 
slope  selected  (vertical  wall  down  to  -5  ft  and  a l-on-1.5  slope)  may  be 
considered  as  a plane  l-on-1.5  side  slope  for  calculation  of  the  reflec- 
tion coefficient. 

53.  In  the  distorted-scale  model,  steepness  of  the  typical  side- 
slope  section  will  be  increased  by  a factor  of  U to  a side  slope  of  1 on 
3/8  as  shown  in  Figure  ll+.  Again  considering  the  distorted-scale  side 
slope  as  a plane  slope,  the  theoretical  reflection  coefficient  is  0.9 
for  small  wave  steepness  in  the  20-sec  and  longer  wave  period  range. 
Thus,  no  modification  in  the  shape  of  the  distorted-scale  side  slope  is 
necessary  for  simulation  of  prototype  reflection  in  the  distorted-scale 
model . 

5^.  The  surface  of  the  typical  side  slope  modeled  in  the  two- 

dimensional  wave  tests  was  smooth.  In  the  harbor,  many  side  slopes  are 

covered  by  a protective  layer  of  rock.  However,  the  size  of  the  rock  on 

the  side  slopes  is  small  relative  to  the  wavelength.  Experimental 
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data  with  various  relative  roughness  values  indicate  that  the  slope 
protection  will  have  little  effect  on  the  reflection  coefficient,  and 
the  slope  protection  was  not  reproduced  in  the  model. 

Wave  Transmission 


55.  Long-period  wave  energy  may  enter  a harbor  either  by  penetra- 
tion through  rock-fill  breakwaters  or  diffraction  through  harbor  en- 
trances. Wave  penetration  or  transmission  through  the  harbor  breakwater 
is  especially  important  in  the  case  of  Los  Angeles  and  Long  Beach  Har- 
bors due  to  the  long  extent  of  the  breakwater.  Energy  transmission 
through  the  breakwater  will  affect  harbor  resonance  and  should  be  simu- 
lated in  the  model.  Two-dimensional  flume  tests  were  conducted  to 


estimate  the  wave  transmission  coefficient  of  a typical  section  of  the 
Long  Beach  breakwater  (illustrated  in  Figure  5 with  the  top  of  the  im- 
permeable section  at  an  elevation  of  -2h  ft).  Typical  prototype  dimen- 
sions (diameter  D)  of  the  Class  A and  Class  B rock  are  and  1.8  ft, 
respectively.  The  breakwater  section  was  modeled  at  a 1:20  scale  with  a 
prototype  water  depth  of  50  ft.  Wave  transmission  tests  were  conducted 
in  a 2-ft-wide  flume  approximately  150  ft  long  with  a flapper-type  wave 
generator.  The  wave  transmission  test  series  was  conducted  in  the  20- 
to  100-sec  wave  period  range.  The  wave  period  range  was  limited  to 
100  sec  by  the  wave  flume  length.  The  wave  transmission  coefficient  C^ 
was  calculated  from 

^t  = h 

1 

where  is  the  transmitted  wave  height.  The  calculated  transmission 

coefficient  for  the  20-  to  100-sec  period  range  is  shown  in  Plate  25  as 
the  curve  for  undistorted-scale  conditions.  The  transmission  coeffi- 
cient varies  from  0.17  for  20  sec  to  0.2h  for  100  sec. 

56.  A wave  transmission  test  series  also  was  conducted  in  the  two- 
dimensional  flume  for  three  distorted-scale  breakwater  sections  with  a 
1:100  vertical  scale  and  a l:i*00  horizontal  scale  as  shown  in  Figure  15. 
The  three  breakwaters  tested  included  a cross  section  composed  entirely 
of  the  larger  (Class  A)  rock  scaled  by  the  vertical  scale  ratio,  a cross 
section  composed  entirely  of  the  smaller  (Class  B)  rock  scaled  by  the 
vertical  scale  ratio,  and  a cross  section  with  Class  A and  Class  B rock 
scaled  by  the  vertical  scale  ratio.  The  calculated  wave  transmission 
coefficient  for  each  of  these  three  cross  sections  is  shown  in  Plate  25- 
As  shown  by  the  wave  transmission  data  in  Plate  25,  the  section  composed 
only  of  the  larger  (Class  A)  rock  has  a higher  transmission  coefficient 
over  the  period  range  than  the  undistorted-scale  cross  section,  and  the 
section  composed  only  of  the  smaller  rock  has  a lower  transmission  coef- 
ficient. However,  the  wave  transmission  coefficient  for  this  breakwater 
section  composed  of  Class  A and  B rock  agrees  closely  with  the  data  for 
the  undistorted-scale  section.  Therefore,  this  breakwater  section  was 
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Figure  15.  Typical  section  of  model  break- 
water; distortion  = ^4;  rock  sizes  scaled 
using  vertical  scale  (l:100) 


selected  for  installation  in  the  model  for  the  Middle  and  Long  Beach 
breakwaters. 

57.  As  shown  in  Figure  5»  the  San  Pedro  Breakwater  has  a solid 
granite  cap  but  does  not  have  an  impervious  core.  The  rock  size  is 
Class  A or  a typical  dimension  of  I4.7  ft.  The  vertical  scale  ratio  of 
1:100  also  was  used  for  selecting  rock  in  the  San  Pedro  breakwater. 
Without  the  impervious  core  or  the  smaller  (Class  B)  rock,  the  San  Pedro 
breakwater  should  have  a larger  transmission  coefficient  than  the  other 
two  breakwaters.  Using  the  following  equation  developed  by  Dr.  G.  H. 
Keulegan  for  wave  transmission  in  the  prototype 


(15) 


where  AL  is  the  length  through  the  breakwater  structure;  the  wave 


transmission  coefficient  for  a 1-ft  wave  height  and  an  average  length  at 
middepth  of  128  ft  is  0.63.  Equation  15  is  valid  only  for  long  waves 
where  the  wave  celerity  may  be  approximated  by  . Equation  15  indi- 

cates the  transmission  coefficient  will  increase  as  the  wave  height  de- 
creases. The  complete  equation  without  the  long  wave  approximation  can 
be  found  in  Reference  26. 

Wave  Absorbers  and  Filters 


58.  Wave  absorber  in  the  Los  Angeles-Long  Beach  Harbors  model  is 
composed  of  2-in. -thick,  2-ft-wide  layers  of  rubberized  hair.  Three 
adjacent  rows  of  rubberized  hair  are  stacked  along  the  perimeter  of  the 
model  for  a total  absorber  width  of  6 ft.  In  a previous  study,  the 
diameters  of  the  hairs  were  found  to  be,  for  the  most  part,  approxi- 
mately 0.00051  ft  and  the  rubberized  hair  material  had  a porosity  of 

0.96.  The  wave  height  transmitted  through  the  filter  may  be  calculated 
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from  the  following  equations  developed  by  Dr.  G.  H.  Keulegan : 


+ B 


(16) 


where 

= initial  wave  height  in  the  filter 
= final  wave  height  before  leaving  the  filter 
AL  = length  of  the  filter 

B = 5.i4  (LT)^^^  (IT) 

for  a kinematic  viscosity  of  0.9  (10~^). 

59-  From  Equations  I6  and  17  and  assuming  a reflection  coeffi- 
cient of  1.0  for  the  model  wall,  the  ratio  of  the  incident  wave  height 
to  the  final  wave  height  for  a 0.05-ft  model  wave  height  transmitted 
through  the  rubberized  hair  to  the  model  wall  and  reflected  back  through 
t le  hair  again  in  model  depths  of  0.5  and  3.0  ft  is: 


Model  Wave  


Period,  sec  Depth  3.0- ft  Depth 

0.6U  0.01  0.06 

1.09  0.01  0.09 

1.56  0.02  0.10 

2.03  0.02  0.11 

2.52  0.02  0.11 

3.7i+  0.02  0.12 

5.01  0.02  0.12 

10.0  0.02  0.12 


Equation  l6  has  been  confirmed  with  experimental  data  in  the  O.TT-  to 
1.91-sec  wave  period  range  but  has  not  been  confirmed  with  experimental 
data  in  the  longer  period  range.  The  preceding  reflection  data  indi- 
cate that  6 ft  of  wave  absorber  along  the  model  perimeter  will  be  suffi- 
cient. Provision  has  been  made  in  the  model  for  an  additional  !+-ft 
thickness  of  wave  absorber  if  required  for  future  wave  tests.  The  wave 
will  also  be  partially  reflected  from  the  front  surface  of  the  rubber- 
ized hair  as  the  wave  first  impinges  on  the  hair,  but  WES  experience 

with  rubberized  hair  wave  absorber  in  a previous  hydraulic  model  in- 
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vestigatio'^  of  oscillation  in  Port  Hueneme  Harbor,  California,  indi- 
cates that  surface  reflections  will  not  adversely  affect  the  generated 
wave  heights. 

60.  In  the  tidal  headbay  area  of  the  model,  the  width  of  the  wave 
absorber  has  been  decreased  to  ft.  The  decrease  in  width  will  allow 
for  tidal  inflow  during  circulation  studies  without  modification  of  the 
absorber.  fhe  remaining  wave  absorber  will  not  adversely  affect  normal 
tidal  circulation. 

61.  Two-dimensional  wave  flume  tests  also  were  conducted  for 
design  of  wave  filters  to  be  placed  in  front  of  the  wave  generator  for 
attenuation  of  waves  reflected  from  offshore  topography  and  the  harbor 
breakwater.  The  wave  filter  test  series  was  conducted  in  the  2-ft 
flume  with  a flapper-type  wave  generator.  The  filters  were  made  using 
sections  of  2-in. -thick  rubberized  hair  retained  in  metal  frames  and 
suspended  on  edge  in  the  flume.  The  frames  were  made  from  angle  iron 
with  a 6-in. -square  wire  mesh  holding  the  hair  in  place.  The  frames 
were  1 ft  apart  in  the  flume  and  tests  were  run  with  2,  h,  and  6 filter 


frame  sections  placed  in  series.  The  wave-height  transmission  coeffi- 
cients for  the  test  configuration  are  shown  in  Figure  l6  for  the  0.5-  to 
10-sec  model  period  range.  The  data  show  the  transmission  coefficients 
are  small  at  0.5  sec  but  rise  rapidly  to  approximately  3 to  ^4  sec  and 
then  continue  with  a gradual  increase  in  the  remaining  period  range. 

The  data  for  two  filter  sections  in  series  are  shown  only  in  the  0.5-  to 
i4-sec  range. 

62.  Arrangement  and  choice  of  the  number  of  filter  sections  to 
be  used  will  be  based  somewhat  on  experience  gained  during  model  opera- 
tion. Initially,  it  is  anticipated  that  in  the  shorter  period  range 
(up  to  approximately  U5- %to  60-sec  prototype),  the  wave  travel  time 
between  the  generator  and  the  breakwater  will  be  sufficient  to  allow 
time  for  any  resonant  oscillations  to  form  in  the  harbor  and  for  measure- 
ments to  be  taken  prior  to  arrival  of  any  re-reflected  waves  from  the 
wave  generator.  Therefore,  wave  filters  may  not  be  needed  in  this 
period  range..  Above  this  period  range,  however,  various  combinations  of 
the  filter  sections  will  be  required. 
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Figure  l6.  Wave  transmission  coefficient  for  rubberized 
hair  filter  sections 
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PAST  III:  MODEL  WAVE  GENERATOR  OPERATION 

AND  PERFORMANCE  REQUIREMENTS 


63.  In  a physical  model  study  of  resonant  response  in  Los  Angeles 
and  Long  Beach  Harbors,  the  wave  generator  should: 

a.  Be  capable  of  generating  waves  with  a prototype  period 
ranging  from  approximately  15  to  6OO  sec. 

Have  a small  variation  in  the  generated  wave  period  and 
height  between  wave  generator  unit  to  unit. 

c_.  Be  capable  of  defining  the  response  over  a narrow  period 
band. 

Be  capable  of  generating  a variable  wave  height  along  the 
wave  front . 

£.  Be  capable  of  generating  a curved  wave  front. 

6U . In  the  period  range  where  the  time  scale  is  approximated  for 
shallow-water  conditions  by  Equation  Tb,  the  time  scale  ratio  is  constant 
at  l:Lo.  Variation  in  the  generated  model  wave  period  of  0.01  sec  would 
correspond  to  a prototype  variation  of  0.1^0  sec.  Results  of  a finite 
element  analysis  of  resonant  response  in  the  Los  Angeles  and  Long  Beach 
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Harbors  ’ indicate  that  the  calculated  maximum  wave-height  amplifica- 
tion may  occur  over  a relatively  narrow  period  band  as  shown  in  Fig- 
ure 17  (adapted  from  Reference  30)  for  prototype  data  gage  5 located  in 
the  southeast  basin  near  Berth  232.  A variation  in  the  generated  model 
wave  period  near  the  maximum  response  may  tend  to  mask  the  actual  maximian 
response  amplitude.  Similarly,  variation  in  the  generated  wave  height 
will  affect  the  level  of  incident  wave  energy  in  the  harbor  and  could 
affect  the  maximum  resonant  amplitude.  Generally,  past  WES  studies  such 
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as  the  hydraulic  model  investigation  of  Port  Hueneme,  California,  in- 
dicate that  the  physical  model  tends  to  respond  over  a broader  period 
band  and  at  a lower  amplitude  where  a sharply  peaked  maximum  resonant 
response  is  calculated  numerically.  The  wave  generator  should,  however, 
have  a sufficiently  small  generated  wave  period  variation  to  adequately 
define  the  response  within  a relatively  narrow  period  band. 

65.  The  necessity  for  variation  of  the  generated  wave  height 
along  the  wave  front  and  a curved  wave  front  has  been  discussed  in  the 
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Figure  17-  Calculated  wave-height  amplification  at  prototype 
wave  gage  location  5 from  numerical  harbor  oscillation  analysis 
(adapted  from  Reference  30) 

wave  refraction  analysis  section  of  Part  II.  Prototype  wave  heights 

over  the  low  frequency  portion  of  the  15-  to  600-sec  period  range  are 

quite  small  and  could  not  be  measured  accurately  in  the  model;  therefore, 

increased  wave  heights  will  be  used  during  the  model  test  series  and  the 

wave  response  data  will  be  normalized.  The  increase  in  wave  height 

should  be  limited  to  the  wave-height  range  in  which  the  increased  height 

will  have  little  effect  on  the  wave  form. 

66.  In  evaluating  the  effect  of  an  artificial  increase  in  wave 
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height,  the  second  order  Stokes  wave  theory,  where  the  ratio  of  water 
depth  to  wavelength  is  1/10  or  larger,  may  be  used;  and  for  a depth  to 
wavelength  ratio  of  less  than  1/10,  cnoidal  theory  may  be  used.  The 
second  order  Stokes  theory  predicts^^  that  for  a finite  amplitude  wave, 
the  amplitude  of  the  crest  is  larger  than  the  trough  amplitude  and  the 
form  of  the  wave  crest  is  more  peaked  and  shorter  than  the  trough.  At  a 
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depth  to  wavelength  ratio  of  1:10  (prototype  period  of  approximately 
UU  sec),  the  ratios  of  the  wave  crest  amplitude  a^  to  the  wave  height 
for  height  ratios  of  0.025h,  0.05h,  O.lOh,  and  0.20h  where  h is  equal 
to  the  average  depth  of  39  ft  are: 


Wave-Height 

Ratio 


Crest  Amplitude  to 
Wave-Height  Ratio 


0.025  0.52 

0.05  r ‘ 0.>3 

o.iO‘  , 0.56 

0.20  0.62 


The  preceding  variation  will  be  the  maximum  in  the  range  h/L'  10  f9r 
the  second  order  theory. 

67.  For  a prototype  period  of  100  sec,  the  h/L  ratio  is  0.17  for 
the  39- ft  average  depth.  The  wave  crest  amplitude  to  wave-height  ratios, 
again  for  wave  heights  of  0.025h,  0.05h,  O.lOh,  and  0.20h  where  h 
equals  the  average  harbor  depth  are: 


Wave-Height 

Ratio 


Crest  Amplitude  to 
Wave-Height  Ratio 


0.025 

0.05 

0.10 

0.20 


0.  56 
0.62 
0.70 
0.78 


At  longer  periods,  the  wave  crest  is  more  sharply  peaked  (smaller  wave 
trough  amplitude)  for  each  wave-height  ratio. 

68.  The  effects  of  finite  amplitude  on  wave  celerity  may  be 
approximated^^  by 


(18) 


for  waves  of  long  length  where  the  curvature  of  the  profile  is  small. 

For  the  depth  to  wave-height  ratio  of  1:10,  the  ratios  of  the  wave  ce- 
lerity to  the  shallow-water  approximation  of  wave  celerity  (Equation  5) 
for  wave  heights  of  0.025h,  0.05h,  O.lOh,  and  0.20h  where  h equals  the 
average  harbor  depth  are: 
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Wave-Height 
Ratio 


Celerity 

Ratio 


0.025 

1.01 

0.05 

1.02 

0.10 

i.oi* 

0.20 

1.07 

The  effect  of  the  finite  amplitude  increases  as  the  wave  amplitude  in- 
creases and  has  more  effect  on  the  surface  wave  form  than  on  wave  ce- 
lerity. In  order  to  minimize  the  effects  of  finite  amplitude,  the  maxi- 
mum wave  heights  for  long-period  waves  in  the  harbor  area  should  be 
limited  to  0.20h  and,  where  possible,  less  than  O.lOh.  For  the  1:100 
vertical  and  average  depth  of  39  ft,  the  maximum  wave  height  in  the 
model  thus  should  be  limited  to  0.078  ft. 

69.  As  a result  of  the  preceding  operational  criteria,  an  electro- 
hydraulic  wave  generator  was  selected  for  the  model  study.  The  wave 
generator  is  composed  of  l^-i  units  with  15-ft  wave  paddles.  The  sections 
may  be  positioned  to  approximate  a curved  wave  front,  as  indicated  in 
Plates  9-11 , and  the  wave  height  may  be  varied  for  each  wave  generator. 

A 15-ft  unit  of  the  wave  generator  with  the  frame,  wave  paddle,  and 
hydraulic  power  supply  is  shown  in  Figure  I8. 

TO.  A mechanical  wave  generator  with  an  electric  motor  and  me- 
chanical linkage  similar  to  those  presently  in  use  at  WES  was  consid- 
ered; however,  the  necessity  for  small  variations  in  the  generated  wave 
period,  difficulty  in  setting  an  extremely  exact  wave  paddle  stroke,  and 
the  requirement  for  model  wave  periods  near  0.5  sec  precluded  their  use. 
Also,  a curved  wave  front  with  wave-height  variation  along  the  wave  front 
would  be  difficult  to  obtain  with  a mechanical  wave  generator.  A pneu- 
matic wave  generator  was  considered  for  use  in  the  model  study  and  sev- 
eral modifications  of  a pilot  pneumatic  wave  generator  were  constructed 
and  evaluated  in  a two-dimensional  wave  flume  ( 10  ft  wide  with  a 3- ft 
water  depth).  During  the  evaluation,  oscillations  of  the  pneumatic 
chamber  produced  a secondary  wave  of  approximately  0.8  sec  on  the 
primary  wave  form  above  approximately  2 sec.  The  height  of  the 
secondary  wave  was  of  the  same  order  of  magnitude  as  the  wave  height 
to  be  used  in  the  model  study.  Several  methods  (principally  additional 
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Figure  l8.  Electrohydraul ic  wave  generator  section  with  frame,  wave 

paddle,  and  hydraulic  supply 

frame-bracing,  baffles,  and  filters)  were  tried  in  an  effort  to  eliminate 
the  secondary  wave  but  were  not  successful.  It  became  apparent  that  a 
substantial  research  and  development  effort  was  required  to  overcome 
these  problems.  Consequently,  the  pneumatic  generator  approach  was  aban- 
doned for  this  project  model  investigation.  The  preceding  evaluation  is 
valid  only  for  the  pilot  pneiimatic  generator  configuration  considered  by 
WES  and  for  the  rather  large  operating  range  required  and  should  not  be 
considered  as  generally  applicable  to  pneumatic  wave  getierators. 

71.  As  a result  of  the  preceding  pneumatic  wave  generator 
analysis,  and  consieiering  the  advances  made  in  electrohydraul ic  systems 
in  recent  years,  the  electrohydraul ic  wave  generator  concept  was 
selected  for  evaluation.  Initially,  a pilot  electrohydraul ic  wave  gen- 
erator was  designed  and  constructed  at  WliS  using  a digitally  controlled 
actuator  (hydraulic  cylinder  used  to  move  the  wave  paddle)  with  a 
position  feedback  control  system.  An  evaluation  of  the  pilot  generator 
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with  a digitally  controlled  actuator  indicated  that  the  electrohydraulic 
wave  generator  was  feasible  but  the  digitally  controlled  actuator  did 
not  provide  sufficient  velocity  response  in  the  high-frequency  range. 

To  improve  the  actuator  response,  the  final  wave  generator  design  in- 
cluded an  electrohydraulic  actuator  with  an  analog  servocontroller  as 
shown  in  Figure  19.  The  main  components  in  each  drive  system  are  the 


Figure  19-  Schematic  diagram  of  electro- 
hydraulic actuator  system 

actuator  with  servovalve  and  linear  voltage  differential  transformer 
(LVDT),  servocontroller,  hydraulic  power  supply,  and  interconnecting 
electrical  cables  and  hoses.  The  electrohydraulic  actuator  systems  were 
supplied  by  a private  contractor  and  the  wave  generator  frames  and  wave 
paddles  were  designed  and  built  at  WES.  Performance  requirements  for 
each  of  the  1^  actuator  systems,  including  the  controllers,  were  that 
each  actuator  system  should; 

a.  Respond  over  a frequency  range  from  0.001  to  2 Hz. 

b.  Have  a displacement  of  jf6  in.  at  0.001  Hz  and  ^,1.3  in.  at 
2 Hz  in  response  to  a sinusoidal  command  input  signal. 
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c^.  Be  capable  of  responding  to  any  command  signal  frequency 
and  amplitude  which  does  not  exceed  the  ^ in.  displace- 
ment limit  or  the  velocity  limit  corresponding  to 
. 3 5 in.  at  2 Hz. 

Maintain  a static  actuator  position  error  of  3 ess  than 

. 06  in. 

e^.  Maintain  a peak-to-peak  displacement  error  within  1% 
relative  to  the  command  signal  amplitude  over  the 
amplitude  range  of  ^1.2  to  _+6  in.  in  the  frequency  range 
from  0.001  to  2.0  Hz. 

Have  a phase  lag  between  the  electrical  input  control 
signal  and  the  mechanical  actuator  displacement  of  not 
more  than  20  deg  at  maximum  displacement  with  a phase 
lag  variation  from  cycle  to  cycle  of  not  more  than  h^I  deg 

Have  a phase  lag  difference  from  actuator  to  actuator 
not  to  exceed  U deg  at  maximum  displacement. 


In  the  preceding  performance  requirements,  no  specific  requirement  is 
given  for  maximum  frequency  variation;  however,  the  actuator  must  follow 
the  command  signal  displacement  within  the  limits  imposed  by  requirement 
^ and  f^.  The  performance  requirements  for  the  wave  generator  were 
selected  so  that  the  wave  generator  can  be  used  in  other  models  or  wave 
flumes  with  similar  operational  requirements.  The  wave  generator  is 
capable  of  generating  wave  heights  (paragraph  73)  in  excess  of  those 
that  will  probably  be  required  for  the  Los  Angeles-Long  Beach  model; 
however,  the  wave-height  capability  may  be  needed  for  future  studies. 

In  the  0.5-  to  1.0-sec  period  range,  the  maximum  wave  height  will  prob- 
ably be  steepness  limited. 

72.  The  wave  height  and  force  required  to  generate  a particular 
wave  for  the  3-ft  water  depth  in  the  Los  Angeles-Long  Beach  Harbors 
model  can  be  estimated  using  small  amplitude  theory  and  neglecting  the 
effects  of  viscosity,  turbulence,  and  leakage  around  the  wave  paddle. 

The  development  of  the  theory  is  described  in  Reference  32  and  is  sum- 
marized briefly  in  the  following  section.  The  velocity  potential  4> 
for  waves  propagating  in  the  +x  direction  is  given  by 


hso 

a 


-b  cosh  (m  z)e 
o o 


im  x 
o 


i L'  b 

n=l 


COS 


(m  z)e 
n 


-m  z 
n 


-iat 
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where  the  center  of  the  coordinate  system  is  at  the  base  of  the  wave 

paddle,  y is  the  positive  upward,  b and  b are  coefficients,  s 

on  ’ 

is  the  amplitude  of  the  displacement  of  the  wave  paddle,  and 


(19a) 


The  coefficient  m is  defined  by 
o 


m K tanh  m h = 
o*^  o 


(20) 


and  m is  defined  by 
n •' 


m g tanh  m h = -a^ 
n*^  n 


(20a) 


The  yelocity  potential  satisfies  the  boundary  conditions  that  the  free 
surface  is  at  a constant  pressure,  the  yelocity  of  the  wave  paddle  and 
the  water  adjacent  to  the  wave  paddle  are  the  same,  and  the  vertical 
velocity  across  the  x-axis  is  zero.  The  velocity  potential  and  its 
first  derivatives  are  finite  in  the  region 


0 < X 


and 


0 < z < h 


and  the  coefficients  b^  and  b^  are  evaluated  through  the  boundary 
condition  that  the  horizontal  wave  paddle  and  water  velocities  be  the 
same. 

73.  The  generated  wave  height  is  given  by 


z = h 

X ->■  00 


or,  having  evaluated  b^  , 


(21) 
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— = 2b  cosh  m h 
so  o 


The  pressure  fluctuation  per  unit  width  along  the  face  of  the  wave 
paddle  may  be  evaluated  from 


P = P [-  ff- 


p = - pgsib^  cosh  m^z  + pgs  ^ b^  cos  (m  z)  e + pg(h  - z)  (23) 

n=l 


The  first  term  in  Equation  23  represents  the  dynamic  pressure  on  the 
wave  paddle  and  is  composed  of  a resistive  part  in  phase  with  the  wave 
paddle  velocity  and  an  inertial  part  in  phase  with  the  wave  paddle 
acceleration.  The  second  term  in  Equation  23  represents  the  hydrostatic 
pressure  on  the  wave  paddle,  and,  in  the  Los  Angeles-Long  Beach  Harbors 
model,  is  canceled  by  the  water  pressure  on  the  rear  of  the  wave  paddle. 

7^*.  The  force  F per  unit  width  is  given  by 


Substituting  from  Equation  23,  Equation  2h  may  be  rewritten  as 


F = -(iFj^  + Fj)e" 


where  F is  the  resistive  force  and  F_  is  the  inertial  force  defined 
n 1 


pgsb^ 

F„  = sinh  m h 

R m o 

o 


“ m h 

Ft  = PSB  E t sin 


The  maximiun  force  F during  a cycle  is 

m 
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(27) 


F 

m 


The  wave  generator  can  be  operated  in  either  a piston  or  fl§,pper  mode 
as  shown  in  Figure  20.  The  solutions  for  and  b^  are  dependent 


► DIRECTION  OF  WAVE-PADDLE  MOVEMENT 

? / rv  / /■  / / / V\  WAVE  PADDLE 

Figure  20.  Wave  generator  piston  and  flapper  modes  of  wave  generation 

on  the  mode  of  operation;  therefore,  the  wave  height  and  maximum  force 
will  vary  with  the  mode  of  operation.  The  maximum  calculated  wave 
height  that  can  be  generated  and  maximum  force  on  the  wavepaddle  are 
shown  in  Figures  21  and  22  over  the  0.5-  to  10.0-sec  period  range.  The 
maximum  wave  height  occurs  near  1.25  sec  for  the  flapper  mode  and  near 
1.75  sec  for  the  piston  mode.  The  decrease  in  amplitude  below  these 
periods  is  due  to  the  decreased  length  of  the  maximum  stroke  as  the 
period  decreases.  As  indicated  by  Figure  21,  the  maximum  calculated 
height  for  the  larger  wave  periods  in  the  flapper  mode  is  approximately 
one-half  the  height  in  the  piston  mode.  Near  0. 5 sec  in  the  short- 
period  range,  the  wave  heights  approach  the  same  value. 

75-  Selection  of  the  wave  generator  mode  of  operation,  either 
flapper  or  piston,  is  somewhat  at  the  discretion  of  the  model  operator, 
provided  a sufficient  wave  height  can  be  obtained  in  the  flapper  mode. 
However,  the  calculated  vertical  profile  over  depth  of  horizontal 
velocity  for  a wave  is  more  nearly  matched  by  the  velocity  regime  gen- 
erated by  the  piston  mode  at  wave  periods  larger  than  approximately 
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MOOEU  WAVE  PERIOD.  SEC 


Figure  21.  Calculated  maximum  wave  height  in  the  piston 

and  flapper  modes 


MODFl  iNAvE  PERIOD.  5tC 


Figure  22.  Calculated  maximum  force  for  each  wave  generator  section 
in  the  piston  and  flapper  modes 
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1.75  sec  in  the  model.  A transition  region  near  1.75  sec  is  suggested 
for  the  change  from  flapper  to  piston  mode.  The  theoretical  wave  veloc- 
ity profile  calculated  from  small  fimplitude  wave  theory  for  1.75  sec  is 
shown  in  Figure  23  for  a 3-ft  water  depth.  The  wave  velocity  is  nor- 
malized by  the  surface  velocity.  Also  shown  in  Figure  23  are  the  veloc- 
ity profiles  for  the  flapper  and  piston  modes.  At  lower  periods,  the 
wave  velocity  profile  tends  to  approach  the  flapper  velocity  profile 
and  at  longer  periods,  the  piston  velocity  profile.  Without  a transi- 
tion point  from  the  flapper  to  piston  mode,  a calibration  of  the  wave 
generator  for  wave  height  versus  stroke  would  be  necessary  for  those 
wave  periods  used  with  both  modes. 

76.  Equations  21  and  27  for  the  wave  height  (except  over  a 

32 

limited  range)  and  wave  force  have  not  been  verified  experimentally. 

32 

Available  data  do  indicate  that  leakage  around  the  wave  generator  may 
significantly  affect  the  generated  wave  height,  and  the  wave  generator 
must  be  calibrated  experimentally  for  wave  height  versus  stroke.  Also, 
initial  data  from  test  runs  with  the  wave  generator  indicate  performance 
requirement  e of  paragraph  71  for  a 1 percent  peak-to-peak  di spilacement 


NORMALIZED  VELOCITY 

Figure  23.  Normalized  theoretical,  flapper  mode, 
and  piston  mode  vertical  wave  velocity  profiles 
for  a 1.75-sec  model  wave  in  a 3-ft  depth 
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error  relative  to  the  command  signal  amplitude  may  be  difficult  to  main- 
tain consistently,  and  this  requirement  may  have  to  be  relaxed  slightly 
during  model  operation. 


Wave  Generator  Frame  Design 


77.  The  vave  generator  support  frame  for  each  unit  is  constructed 
of  structural  grade  steel  with  welded  joints.  The  frame  is  10  ft  long 
by  6 ft  wide  by  3.7  ft  high.  Front  and  end  elevations  of  the  frame 
with  the  wave  paddle  are  shown  in  Figure  2h.  The  frame  is  mounted  on 
swivel  casters  and,  when  in  position,  can  be  leveled  using  screw-type 
adjustable  legs.  The  frame  includes  a drip  pan  to  prevent  contamination 
of  water  in  the  model  in  the  event  of  a hydraulic  fluid  leak. 

78.  The  wave  paddle  is  driven  by  two  2-in.  drive  rods  at  points 
ft  from  the  ends  of  the  wave  paddle.  The  drive  rods  are  chrome- 

coated  stainless  steel  and  each  rod  is  supported  by  two  bronze  bushings 
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Figure  2k.  Sectional  views  of  the  wave  generator  frame 
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mounted  on  an  8- in.  channel  section.  The  bottom  of  the  wave  paddle  is 
supported  by  two  auxiliary  drive  rods  at  points  14.0^4  ft  from  the  ends 
of  tVie  wave  paddle.  The  force  from  the  actuator  is  transmitted  to  the 
drive  rods  by  an  aluminum  crossarm  assembly  shown  in  Figure  2h . In  the 
piston  mode,  the  auxiliary  drive  rod  is  connected  to  the  crossarm  as- 
sembly. In  the  flapper  mode,  the  auxiliary  drive  rod  is  connected  to 
the  support  frame  and  the  wave  paddle  rotates  about  the  auxiliary  drive 
rod  pinned  connection  at  the  wave  paddle. 

T9-  The  wave  paddle  is  15  ft  long  and  3-75  ft  high,  and  is  made 
up  of  2-  by  5-in.  rectangular  aluminum  tubes.  The  tubes  are  open  on 
either  end  to  reduce  the  buoyant  uplift  force  on  the  wave  paddle.  The 
actuator  pushes  against  a reaction  block  welded  to  the  support  frame 
and  the  aluminum  crossarm  assembly  to  drive  the  wave  paddle. 

80.  The  wave  generator  frames  were  initially  painted  with  a lead 
primer  coat  and  a latex  enamel  surface  coat;  however,  they  have  recently 
been  repainted  with  a coal  tar  epoxy  on  the  basis  of  recommendation* 
from  the  U.  S.  Army  Construction  Engineering  Research  Laboratory  (CERL) 
to  provide  additional  corrosion  protection. 


• Telephone  communical Ion  between  Dr.  R.  W.  Whalin  of  WES  and  Mr.  Bob 
Quatroon  of  CERL. 
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PART  IV:  MODEL  DATA  ACQUISITION  AND  CONTROL  SYSTEM 


System  Requirements 

81.  The  extreme  complexity  and  large  size  of  the  Los  Angeles- 
Long  Beach  Harbors  model  as  well  as  the  anticipated  volume  of  model 
data  to  be  collected  have  forced  abandonment  of  conventional  techniques 
of  model  data  acquisition  and  control  and  model  data  analyses.  For  ex- 
ample, manual  procedures  for  collecting  and  analyzing  this  model's  wave 
data  woiild  require  approximately  seven  man-years  to  reduce  and  analyze 
the  anticipated  volume  of  data  without  considering  the  efforts  to  man- 
ually collect  these  data.  To  meet  the  demands  of  this  model,  an  auto- 
mated data  acquisition  and  control  system  (ADACS)  with  supporting  soft- 
ware for  model  control,  data  acquisition,  and  analyses  has  been  designed, 
implemented,  and  checked  out. 

82.  In  the  model,  the  variables  to  be  monitored  during  harbor 
oscillation  studies  are  water-surface  elevation  changes  and  water  ve- 
locity. Based  on  past  experience,  prototype  wave  data,  and  this  model's 
scaling  relations,  the  range  of  model  wave  periods  is  expected  to  be 
0.50  to  10  sec  with  the  minimum  period  of  0.50  sec  posing  the  most 
stringent  constraints  on  ADACS  design.  The  expected  maximum  wave  height 
to  be  monitored  in  the  model  is  about  k.O  in.  (0.333  ft).  Wave  period 
resolution  of  0.01  sec  and  wave-height  accuracy  of  ±0.001  ft  are  desired 
capabilities  of  ADACS  with  ±0.001  ft  accuracy  being  the  most  restrictive 
requirement  for  system  design.  For  design  purposes,  the  upper  limit  of 
the  velocity  variable  is  less  critical  than  the  lower  or  threshold  value. 
Model  velocities  as  low  as  0.05  fps  may  be  of  interest  in  this  model. 

The  desired  threshold  velocity  is  0.01  fps  which  is  a stringent  require- 
ment for  existing  model  velocity  sensors  and  is  very  restrictive  in 
system  design.  In  addition  to  meeting  the  requirements  of  the  above 
variables,  automated  controls  for  electrohydraulic  actuators  (wave 
generator  units)  and  calibration  of  the  wave  gages  are  required  capa- 
bilities of  ADACS. 

83.  The  following  characteristics  and  requirements  were  put 


62 


forth  as  design  criteria  which  must  be  met  by  an  ABACS  for  use  in  the 
Los  Angeles-Long  Beach  hydraulic  model  study. 

a.  Model  wave  characteristics. 

(1)  Wave  frequencies  as  high  as  2 Hz  (wave  period  range 
of  0.50  to  25.00  sec). 

(2)  Wave  heights  in  operating  range  of  0.01  to  0.5  ft. 

(3)  Wave-height  accuracy  of  ±0.001  ft. 

b.  Wave  Gages. 

(1)  Variable  number  of  wave  gages  not  to  exceed  50. 

(2)  Water-surface-piercing,  parallel  wire,  resistance- 
type  gages. 

£.  Sampling  techniques. 

(1)  Data  collected  over  at  least  20  wave  periods. 

(2)  Sampling  frequency  variable  and  high  enough  to  define 
the  first  three  harmonics  of  the  2-Hz  wave  frequency 
(minimum  sampling  rate  of  60  samples/cycle). 

(3)  Minimum  time  delay  (not  to  exceed  6 msec  between 
sampling  digitally  the  first  and  fiftieth  wave  gage 
during  any  one  scan  of  the  50  gages,  and  this  time 
delay  should  be  constant  and  independent  of  the 
sampling  frequency. 

(^)  Time  interval  between  scans  of  all  gages  should  be 
controlled  to  within  a few  microseconds. 

d.  Recording  modes. 

(1)  Digital  recording  of  data  from  all  50  channels  in 
binary  code  with  provisions  for  BCD  recording  of 
specific  information  regarding  test  identification 
and  data  analysis. 

(2)  Digital  data  recorded  on  9-track  magnetic  tape  with 
IBM  compatible  record  format. 

(3)  Continuous  analog  recording  of  50  channels. 

(i*)  Time  correlation  of  digital  and  analog  recording 
modes. 

£.  Calibration  of  wave  gages. 

(1)  Efficient  and  accurate  means  of  calibrating  the  wave 
gages  before  a series  of  wave  tests. 

(2)  Recording  of  calibration  data  in  digital  and  analog 
modes. 
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Wave  generator  controls. 

(1)  Provide  a sinusoidal  analog  voltage  whose  amplitude 
can  vary  within  ±10  volts  full  scale  (F.S. ) at 
0.05  percent  F.S.  increments  and  whose  period  (fre- 
quency) can  change  over  0.50  to  25.00  sec  at  0.01- 
sec  increments. 

(2)  Capable  of  varying  the  number  of  cycles  or  total  time 
of  control  signal  for  generator  units. 

Mobility  of  system. 

(1)  Entire  data  acquisition  system  to  be  sufficiently 
mobile  to  service  various  hydraulic  models  from  dif- 
ferent locations  with  a minimum  effort  and  cost  in 
moving  and  installing  the  system. 

(2)  Compact  interfacing  and  panels  for  wires  to  wave 
gages  in  model,  thus  allowing  efficient  mobility  of 
wave  gages  for  expedient  relocation  within  the  model. 

Environment  for  systeni. 

(1)  ABACS  is  housed  in  a controlled  environment  whose 
temperature  ranges  from  65  to  75°F  and  humidity 
ranges  from  1^0  to  55  percent. 

(2)  Emergency  notification  of  appropriate  personnel  if 
temperature  and  hvutiidity  change  beyond  specified 
range . 


Svstem  Configuration 


8U.  The  system  configuration  (Figure  25)  of  ABACS  consists  of 
the  following  subsystems ; 

Bigital  data  recording  and  controls, 
b.  Analog  recorders  and  channel  selection  circuits. 
c_.  Wave  sensors  and  interfacing  equipment, 
d.  Wave  generator  units  and  control  equipment. 

85.  The  digital  data  recording  and  control  subsystem  is  basically 
a minicomputer  (l  psec  memory  cycle  time) ‘with  the  following  charac- 
teristics and  peripheral  devices: 

a.  32K,  16  bit  words  of  core  memory. 

Analog  to  digital  pack  featuring  6h  analog  inputs 
(±10  volts  F.3.),  U5  kHz  multiplexer,  and  12  bit  (in- 
cluding sign)  analog  to  digital  (A/B)  converter. 
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Figure  25.  Automated  data  acquisition  and  control  system  (ADACS)  for 

hydraulic  wave  models 


c_.  Interval  timer  with  a 1 ysec  counter. 

;d.  A long  line  adapter  and  a universal  controller. 

£.  IVo  direct  memory  access  channels. 

Moving  head  disc  with  one  removable  platter  and  one  non- 
removable platter  (l.l  million  word  storage  capacity, 

100  kc  transfer  rate). 

£.  Magnetic  tape  controller  with  two  9-track  magnetic  tape 
units  (25  in. /sec,  8OO  bpi,  10  kc  transfer  rate). 

ii.  Teletype  unit  with  keyboard/printer  and  10  cps  paper 
tape  reader/punch. 

_i.  96  sense  and  control  lines. 

J_.  Matrix  electrostatic  printer/plotter  (print:  162O  lines/ 

min,  132  char. /line,  and  128  char,  set;  plot:  i*  in. /sec 

and  100  styli/in.). 

86.  The  analog  recording  subsystem  acts  as  a backup  for  ABACS 
and  a visual  display  for  operator  Inspection  of  analog  signals  from  wave 
sensors.  This  subsystem  has  manual  or  automated  selection  and  control 
of  five  12-channel  oscillographs. 

87 . The  wave  sensor  subsystem  includes  the  following  major  com- 
ponents: 

a.  50  wave-height  sensors  and  stands. 

b.  50  channels  of  signal  conditioning  equipment. 

£.  Power  supplies. 

Manual  and  automatic  calibration  equipment. 

88.  The  last  subsystem,  wave  generator  and  controls,  is  composed 
of  the  following  major  components. 

a.  ih  wave  generator  frames,  crossarm  assemblies  with  drive 
rods,  and  wave  paddles. 

b.  lU  electrohydraulic  actuators  with  linear  voltage  diffei- 
ential  transformers,  hydraulic  power  supplies,  and  hy- 
draulic controls. 

£.  1*4  servovalves  and  servocontrollers . 

d.  Interface  from  servocontrollers  to  ABACS. 

Each  unit  of  the  wave  generator  (refer  to  Part  III)  utilizes  the  elec- 
trohydraulic actuator  for  driving  the  wave  paddles  through  the  crossarm 
assembly  and  drive  rods.  The  electrohydraulic  actuator  is  controlled 
by  the  servovalve  and  servocontroller  which  accepts  ABACS  command 
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sit^iidls  controlling  wave  heights  and  jieriods  (rel’er  to  paragraiiiis  9i-95}~ 
Except  for  the  wave  sensors  and  the  wave  generator  units,  all  components 
of  ALAC:'  are  compactly  iioused  in  a trailer  (Figure  26)  to  provide  easy 
control  of  tiie  environment  of  the  system. 


Figure  26.  Automated  data  acquisition  and  control  system 


Wave  Sensors 


89.  'fhe  data  acquired  from  wave  models  are  the  water-surface 
variations  about  a reference  water  level . This  information  is  collected 
at  selected  geographic  locations  within  the  model  for  specified  wave 
conditions  at  the  wave  generator.  Wave  sensors  are  used  to  obtain  this 
information  at  selected  locations  in  the  model.  Water-surface-pi ercing, 
parallel-rod  wave  gages  (Figure  27)  are  the  wave  sensors  used  in  the 
Los  Angeles-Long  Beach  Harbors  model.  This  type  of  sensor  has  been 
used  on  j^ast  wave  models  with  much  success.  However,  for  the  present 
case,  this  sensor  required  some  modifications  to  the  electronics  (to 
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Figure  27.  Parallel-rod  wave  sensor 


increase  accuracy  and  stability)  and  to  the  physical  configuration  (for 
automated  calibration). 

90.  Each  water-surface-piercing,  parallel-rod  wave  sensor  is 
connected  to  a Wheatstone  bridge  (Figure  28).  The  transducer  measures 
the  conductance  of  the  water  between  two  parallel  rods  moutited  verti- 
cally. The  conductance  is  directly  proportional  to  the  depth  of  sub- 
mergence of  the  two  rods  in  the  water.  The  output  of  each  wave  sensor 
is  routed  through  shielded  cables  to  its  signal  conditioning  equipment 
where  it  is  processed  for  recording.  The  output  of  the  signal  condition- 
ing equipment  is  connected  through  shielded  cables  to  analog  oscillo- 
graphs where  an  analog  tim«^  history  is  recorded  and  to  the  analog  multi- 
plexer of  the  digital  recoi'ding  subsystem  where  it  is  digitized  and 
recorded  in  a binary  format  on  magnetic  tape  and/or  disc.  The  signal 
conditioning  equipment  (Figure  29)  consists  of  a carrier  amplifier  and 
various  power  supj'lies.  This  system  can  detect  changes  in  water-surface 
elevations  to  an  accuracy  of  ±0.001  ft.  To  obtain  this  accuracy, 
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P’igure  29-  Signal  conditioning  equipment 
for  wave  rod  amplifiers 


ultrastable  power  supplies  and  better-than-average  signal  to  noise 
ratios  are  necessary.  The  carrier  source  for  the  wave  sensor  bridge 
maintains  a variation  of  less  than  0.025  percent  for  a 10  percent  power- 
line variation.  The  noise  in  the  system  is  less  than  ±5  mv  for  ±10  v 
full  range.  Thus,  the  signal  to  noise  ratio  is  -66  db.  To  maintain 
this  signal  to  noise  ratio,  it  is  necessary  to  use  shielded  cables 
wherever  possible,  to  discriminate  selections  of  grounds,  to  use  high- 
quality  components,  and  to  have  a good  maintenance  schedule. 

91.  To  convert  the  water-elevation  data  in  millivolts  to  water- 
surface  elevations  in  feet,  each  wave  sensor  must  be  calibrated.  The 
capability  of  automatically  calibrating  the  wave  sensors  (maximum  of 
25  rods  simultaneously)  prior  to  collecting  data  also  is  provided  by 
ABACS.  In  order  to  calibrate  each  set  of  parallel  rods,  the  voltage 
from  the  signal  conditioning  equipment  is  monitored  and  recorded  as  the 
parallel  rods  are  moved  vertically  a known  distance  into  or  out  of  the 
water.  A precision,  linear-position  potentiometer  is  located  on  the 
wave  sensor  stand  and  is  coupled  directly  to  the  parallel  rods  by  a 
gear-train  driven  electric  motor.  By  moving  the  coupled  wave  sensor  and 
potentiometer  wiper  vertically  with  the  electric  motor  and  by  monitoring 
the  output  voltage  from  the  potentiometer,  the  wave  sensor  can  be  moved 
a precise  distance.  The  electric  motor  for  each  wave  sensor  is  con- 
trolled by  a control/sense  line  and  a relay  contact.  The  minicomputer 
controls  the  vertical  movement  of  each  wave  s<»nsor  by  actuating  its 
associated  control/sense  line.  The  central  processing  unit  (cpu)  acts 
as  a voltage  comparator  by  monitoring  the  potentiometer  voltage  and 
comparing  it  with  a reference  voltage  which  is  determined  from  desired 
displacement  and  potentiometer  calibration.  When  the  voltage  compari- 
son is  satisfied,  the  control/sense  line  is  reactivated,  the  electric 
motor  stops,  and  voltage  samples  from  the  rods  and  potentiometers  are 
acquired. 

92.  By  systematically  moving  each  wave  sensor  through  11  quasi- 
equally  spaced  locations  (Figure  30)  over  the  range  of  the  rod  length 
used,  voltage  versus  known  displacements  are  obtained  from  which  a 
calibration  curve  for  each  sensor  can  be  calculated  and  recorded  on 
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Figure  30.  Schematic  of  wave  rod  calibration 
procedure 

magnetic  tape  or  disc.  After  collecting  the  calibration  data,  the  mini- 
computer analyzes  these  data  by  least-squares  fitting  a set  of  curves 
(linear,  quadratic,  or  spline)  to  the  data,  determining  the  best  order 
of  fit,  and  comparing  the  maximum  deviation  of  the  best  fit  with  a pre- 
viously acceptable  value  for  this  maximum  deviation.  If  the  fitted 
curves  are  not  acceptable,  the  minicomputer  flags  that  channel  in  the 
calibration  record  on  the  magnetic  tape  or  disc  for  further  analyses. 

Any  malfunctioning  sensors  are  listed  on  the  teletype  for  the  operator 
to  determine  the  required  action  (i.e.,  accept  present  calibration, 
clean  bad  rods,  recalibrate,  etc.).  Having  completed  the  calibration 
mode,  the  original  calibration  data  for  each  set  of  parallel  rods  and 
potentiometer  as  well  as  the  calibration  coefficients  are  written  into 
a file  on  magnetic  tape  or  disc.  Analytical  cons iderations  of  calibra- 
tion curve  fitting,  accuracy  of  reference  potentiometers,  and  repeatabil- 
ity test  of  calibration  procedures  have  demonstrated  that  an  accuracy 
of  ±0.001  ft  in  wave  height  is  obtainable. 


Data  Acquisition 


93.  During  the  acquisition  mode,  wave  data  for  a specified  wave 
condition  at  the  wave  generator  are  collected  from  a maximum  of  50  wave 
sensors,  recorded  on  analog  strip  charts,  digitized,  and  recorded  on 
magnetic  tape  or  disc  for  further  analyses.  In  addition,  actuator  dis- 
placement data  from  the  1^  electrohydraulic  wave  generator  units  are  ac- 
quired. The  sampling  scheme  is  quite  flexible  and  can  be  tailored  for 
different  applications  with  maximum  thru-put  rates  theoretically  limited 
by  the  multiplexer  rate  (U5  kHz)  and  allocatable  buffer  size.  However, 
for  specific  types  of  tests,  the  execution  time  of  applications  software 
and  the  specified  number  of  discrete  samples  per  wave  period  reduce  the 
maximum  thru-put  rate.  The  present  sampling  scheme  used  for  each  wave 
sensor  is  60  discrete  voltage  samples  equally  spaced  over  each  wave 
period  for  a predetermined  number  of  periods.  The  minicomputer  calculates 
from  input  parameters  the  required  timing  interval  between  multiplexer 
scans  to  provide  the  correct  sampling  rate  and  initializes  counters  for 
determining  completion  of  wave  tests.  In  addition,  it  controls  the 
starting  of  the  wave  generator  units,  generates  and  updates  the  control 
signal  to  the  wave  generator  units,  lags  the  beginning  of  data  acquisi- 
tion by  a specified  number  of  wave  periods  after  starting  the  generator, 
and  provides  timing  pulses  for  synchronizing  and  controlling  the  analog 
recorders.  At  completion  of  the  acquisition  mode,  the  calibration  data, 
wave  data,  and  actuator  displacement  data  (paragraph  95)  have  been  re- 
corded in  both  analog  and  digital  form.  These  data  with  a header  for 
test  identification  and  pertinent  parameters  are  available  in  binary  form 
from  disc  for  direct  analyses  by  the  minicomputer  or  on  magnetic  tape  in 
a format  suitable  for  backup  analyses  on  the  WES  Honeywell  0-635  computer. 

Wave  Generator  Controls 


9^.  Through  the  digital  to  analog  converter,  ADACS  provides  to 
the  wave  generator  an  electrical  command  signal  that  drives  the  elec- 
trohydraulic actuator  at  a given  frequency  and  amplitude  of 
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displacement.  The  actuator  drives  the  wave  paddle  which  generates  a 
water  wave  that  has  the  siime  frequency  as  the  actuator  and  a wave  height 
proportional  to  the  actuator's  displacement.  The  control  signal  is  a 
d-c  analog  voltage  that  can  vary  between  110  v.  The  functional  form 
of  the  control  signal  can  be  varied  to  generate  the  required  wave  regime 
such  as  monochromatic,  spectral,  or  any  time-dependent  wave  form.  The 
only  restriction  on  the  form  of  the  control  signal  is  the  maximum  veloc- 
ity limitation  of  the  electrohydraul ic  actuator  system.  For  initial 
wave  tests,  a sinusoidal  or  monochromatic  wave  form  will  be  used,  and 
within  the  desired  range  of  wave  frequencies  and  heights,  no  limitations 
on  actuator  response  exist. 

95.  The  desired  analog  command  signal  from  ABACS  is  applied  to 
the  wave  generator  control  unit  (Figure  3l)  which  provides  lU  parallel 


«AVE  SCALING  SERVOCONTROLLER 

GENERATOR  amplifier  SUMMING  POKER 

CONTROL  UNIT  SPAN  CONTROL  AMPLIFIER  AMPLIFIER 


Figure  31.  Schematic  of  wave  generator  controls 


inverted  output  signals.  These  output  signals  are  adjusted  for  selected 
quiescent  and  amplification  levels  and  then  are  applied  to  the  input  of 
their  respective  scaling  amplifier  of  the  servocontroller . It  is  in 
this  section  of  the  servocontrollt  th  the  span  control  potentiometer 
that  the  system  is  calibrated  for  a +6.000-in.  stroke  for  a +10.000-V 
command  signal.  This  calibrated  command  signal  along  with  the  actuator 
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position  I'eedbfu'k  si*?nul  is  applied  to  the  slimming  junction  of  the  servo 
controller.  The  algebraic  sum  of  these  two  signals  is  then  applied  to  a 
power  amplifier  which  is  suitable  for  driving  the  coil  of  the  servovalve 
The  servovalve  controls  the  flow  of  hydraulic  fluid  to  the  cylinders  in 
the  actuator;  thus,  the  displacement  of  the  piston  is  proportional  to 
the  electrical  command  signal.  The  actuator  position  feedback  trans- 
ducer, a linear  variable  differential  transformer  (LVDT),  provides  a 
modulated  2-kHz  output  which  is  proportional  to  the  actuator  position. 
This  feedback  signal  is  demodulated  and  then  calibrated  for  a scale 
factor  of  approximately  1.66  v/in.  Therefore,  the  feedback  signal  along 
with  the  loop  gain  is  used  to  stabilize  the  closed  loop  electrohydraulic 
system,  thus  ensuring  that  the  actuator  tracks  the  command  input  within 
the  desired  accuracy. 

96.  In  addition,  the  feedback  signal  is  applied  to  the  analog 
digital  converter  (ADC)  of  the  ADACS  where  it  is  used  to  compute  the 
response  of  the  electrohydraulic  system  with  respect  to  the  command  sig- 
nal. From  the  electrohydraulic  actuator  response  (stroke)  and  model 
calibration  data  of  wave  height  versus  actuator  stroke,  the  wave  height 
of  the  generated  wave  for  a specified  frequency  can  be  determined. 

Data  Analysis 

97-  Data  analysis  may  be  performed  by  either  the  minicomputer  or 
by  the  Honeywell  G-635  of  the  Automatic  Data  Processing  Center  at  WES. 
Schematically,  the  procedure  for  wave  data  analyses  is  as  follows: 

a.  Program  initialization. 

(1)  Input  test  parameters  and  option  flags. 

(2)  Read  and  decode  data  files  on  magnetic  tape  or  disc. 

(3)  Demultiplex  data  files  and  scale  data. 

b.  Wave  record  analyses. 

(1)  H ± Ojj  and  T + . 

“rms  ■ 

(3)  H where  x is  a specified  percent  of  the  highest 
wave  heights,  normally  x = 0.333;  thus,  H , is 
significant  wave  height.  ^ 
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(U)  Option  to  plot  wave  heights  versus  time. 
c_.  Fourier  analyses. 

(1)  Autospectrura. 

(a)  Amplitude- frequency. 

(b)  Energy. 

(c)  Spectral  smoothing. 

(d)  Plots  of  above  parameters. 

(2)  Cross  spectrum. 

(a)  Energy. 

(b)  Coherency. 

(c)  Phase. 

(d)  Plots  of  above  parameters, 
d.  Least-squares  harmonic  analysis. 

(1)  Amplitude  and  phase  for  specified  wave  period. 

(2)  Relative  phases  and  amplification  factors  between 

gages. 

(3)  Analysis  of  residual  variance. 

(i*)  Graphic  output  of  above  results. 

98.  The  analysis  of  actuator  displacement  data  is  basically  the 
above  least-squares  harmonic  analysis  from  which  can  be  determined  the 
actuator  response.  The  results  of  the  data  analyses  can  be  plotted  by 
electrostatic  printer/plotter  on  minicomputer  or  by  pen  plotter  or  CRT 
plotter  (microfilm  and  hard  copy)  on  the  Honeywell  G-635  system.  For 
each  model  test,  original  data  and  analyzed  results  of  all  tests  are 
permanently  stored  on  magnetic  tape  and  disc  files  for  future  reference 
and  data  table  generation  for  reports,  etc. 

99.  With  the  analyses  of  wave  records  from  various  harbor  loca- 
tions, amplification  factors  or  harbor  responses  for  various  input  wave 
conditions  are  estimated.  These  results  provide  the  hydraulic  engineers 
with  the  basic  data  whereby  the  effects  of  various  proposed  expansions 
and  modifications  to  existing  harbors  can  be  evaluated. 
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PART  V:  CONCLUSIONS 


100.  Based  upon  the  design  considerations  discussed  in  this  re- 
port, it  is  concluded  that  a 1 ; 100-vertical- , 1 ; UOO-horizontal-scale 
hydraulic  model  of  Los  Angeles  and  Long  Beach  Harbors  will  accurately 
simulate  prototype  wave  conditions  so  that  effects  of  proposed  harbor 
modifications  may  be  accurately  evaluated.  All  predominant  parameters 
(rel'ract  ion,  slioaling,  diffraction,  viscous  friction,  reflection,  trans- 
r.issiori,  l.arbor  oscillation,  incident  waves,  etc.)  have  been  adequately 
re[ reduced  (either  naturally  or  artificially  by  adjustment  of  the  model 
and  its  ajqiurtennnces ) . Tests  of  existing  conditions  should  be  con- 
ducted initially  to  obtain  a base  for  comparison.  The  model  can  then 

be  used  in  a predictive  mode  to  evaluate  various  proposed  modifications, 
thereby  optimizing  the  harbors  design  from  both  a cost  and  a functional 
standpoint . 

101.  Potential  future  uses  of  the  model  not  previously  discussed, 
include  studies  of: 

a.  Dispersion  of  effluent  from  specific  outfalls  in  the 
harbor  complex. 

Tsunami  waves. 

£.  Effects  of  various  coastal  structures  from  Point  Fermin 
to  Huntington  Beach. 

Effects  of  artificial  offshore  islands  or  landfill  sea- 
ward of  the  existing  harbor  breakwater  on  tidal  circula- 
tion and  wave  conditions. 

£.  Tidal  circulation  in  areas  such  as  Huntington  Beach, 
Bolsa  Chica  Wildlife  Reserve,  Alamitos  Bay,  etc. 

£.  Effect  of  various  future  harbor  modifications  on  tidal 
and  wave  phenomena. 

The  physical  environirient  in  the  harbor  and  surrounding 
area. 

]i.  Wave  and  tidal  circulation  in  the  proposed  Downtown 
Marina  and  Queensway  Marina. 
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Table  1 

Comparison  of  Ten  Model  and  Prototype  Orthogonal 
Azimuths  for  a Wave  Period  of  20  sec 


Prototype  Azimuth,  depirees 

2ll*.T 

199. 

197.1 

187.5 

182.8 
178.  i* 

166.8 
188.7 
191.  i* 

216.0 


Model  Azimuth,  degrees 

208.0 

193.0 

192.7 

188.7 
181*.2 

179.8 

167.8 
191.3 
185.2 

225.1 


Table  2 

Comparison  of  Thirteen  Model  and  Prototype  Orthogonal 
Azimuths  for  a Wave  Period  of  75  sec 


Prototype  Azimuth,  degrees 


Model  Azimuth,  decrees 


227.3 

225.9 
217. i* 
209.0 

191.2 
186.7 

178.2 
17i+.6 

155.9 

151.5 

Ul.O 

181. 5 

183.5 


229.7 

226.8 
222.  3 

209.7 

191.3 

181.1 

180.9 

175.7 

153.3 

116.9 

152. 1 
182.  3 

183.1 
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Averaged  value. 


Table  i* 

Wave  Reflection  Cbefficlents  for  the  Distorted-Scale  Middle 
and  Lorip:  Beach  Breakwaters 


itype  Wave 
od,  sec 

Prototype  Wave 
Height,  ft 

Reflection  Coefficient 

20 

0.91 

0.61* 

2.15 

0.50 

30 

1.26 

0,1*1* 

3.01+ 

0.1*1 

1+0 

1.15 

0,1*3 

3.89 

0.1*5 

50 

1.00 

0.1*8 

1+.11+ 

0.51 

60 

1.02 

0.52 

3.75 

0.1*9 

70 

1.03 

0.39 

3.93 

0.1*6 

80 

1.31 

0.1*6 

3.1*6 

0.52 

90 

1.21* 

0.1*8 

l*.6l 

0.  38 

100 

1.55 

0.  39 

3.81 

0. 1*2 

REFRACTION  DIAGRAM  FOR 
30-SECOND  WAVE 
FROM  SOUTH 

GRID  1 


REFRACTION  DIAGRAM  FOR 
60-SECOND  WAVE 
FROM  SOUTH 

GRID  1 


REFRACTION  DIAGRAM  FOR 
30-SEC  WAVES  FROM  SOUTH 
GRID  2 


PLATE  7 


REFRACTION  DIAGRAM  FOR 
360-SEC  WAVES  FROM  SOUTH 
GRID  2 
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BREAKWATER  TRANSMISSION 


APPENDIX  A:  NOTATION 


a Wave  amplitude 

a^  Wave  crest  amplitude 

b ,b  Coefficients  in  velocity  potential  fjinction 
on 

C Wave  celerity 

C Wave  reflection  coefficient 

r 

Wave  transmission  coefficient 
D Typical  rock  dimension 

e Base  of  natural  logarithms 

F Force 

Fj  Inertial  force 

F Maximum  force 

m 

F Resistive  force 

n 

g Acceleration  due  to  gravity 

h Water  depth 

h^  Water  depth  at  deep  end  of  bottom  elevation  change 

h^  Water  depth  at  shallow  end  of  bottom  elevation  change 

H Wave  height 

H Average  wave  height 

Wave  height  at  antinode 
Initial  wave  height 
Wave  height  at  node 
Reflected  wave  height 
Wave  height  at  time  t 
Wave  height  for  closed  harbor 
Transmitted  wave  height 
^RMS  Root-mean-square  wave  height 

Hq  Initial  wave  height  in  wave  absorber 

Final  wave  height  in  wave  absorber  before  leaving  the  absorber 
i Imaginary  number 

k Wave  number 

K Refraction  coefficient 

r 


A1 


Shoaling  coefficient 
L Wavelength 

L Length  scale  ratio 

2 

m Solution  of  mg  tanh  m h = o 

o ° ° 2 

m Solution  of  mg  tanh  m h = -o 

n n°  n 

n Model  distortion 

p Pressure 

R Wave-height  amplification  factor 

s Wave  generator  stroke  amplitude 

t Time 

T Wave  period 

T Average  wave  period 

T Time  scale  ratio 

r 

x,y,z  Coordinate  system  with  x,y  in  the  horizontal  plane  and 
z vertical 

a Dimensionless  parameter 

B Variable 

AL  Width  of  wave  absorber 

n Local  water-surface  elevation 

V Kinematic  viscosity 

p Water  density 

a Angular  wave  frequency 

o„  Standard  deviation  of  the  wave  height 

n 

Orj,  Standard  deviation  of  the  wave  period 

<t>  Velocity  potential 

Subscripts : 

m Model  value 

p Prototype 

r Scale  ratio  of  model  to  prototype 
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